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I. INTRODUCTION 


Numerous empirical equations (see Table I, of which the major por- 
tion is a compilation made by Herrman and Jones!. The geometric values 
listed in the table have been converted to the metric system (cgs) 
system of units) are available for predicting depth of penetration and 
residual velocities for various projectile-target combinations. 
(Projectile and fragment may be used interchangeably through this report.) 
Unfortunately, since the equations are empirical, each equation is 
applicable only in the region of the data used to generate the empirical 
constants. Any extrapolation outside the range of the data is question- 
able. 


Extended use has been made at the BRL of the equations developed 
under the code name "Project Thox''*, The equation for residual velo- 
city for perpendicular impact (0° obliquity) is: 

a boc ,d 

v = Vi 210 (XA) mt, Vs (1) 
The values for the empirical exponents a, b, c and d determined by a 
least squares fit to experimental data are tabulated in Table II for 
various target materials along with the range of plate thickness X,, 
striking velocity V, and the area A and mass of the projectile. 
Since the Thor equation is so widely used at the BRL, in order for any 
model predicting residual velocity to qualify as a replacement to this 
equation, it should be more accurate in its prediction or should exhibit 
other qualities which render it more useful - for example, allow extra- 
polation with greater confidence than a purely empirical model. 


It is desiruble to develop a theoretical model which can be used in 
general to predict terminal ballistics. Needless to say, projectile- 
target interactions are complex. Although kinetic energy projectiles 
have been used for several centuries, no single predictive model has 
been found to be applicable for all test conditions. 


The model to be discussed in this report is completely general in 
that its use decs not require a data base to generate new empirical 
constants. The wnodel is a modification of an analytic approach for 
predicting residual velocity of penetrators impacting targets at 0° 


‘Herrman and Jones, "Correlation of Hyperveloctty Impact Data," 


Proceedings of the Fifth Sympostum on Hyperveloeity Impact, Vol. 1, 
Part 2, April 1962, 


*Projeat Thor, "The Reatstance of Various Metallic Matertale to Perfora-~ 
tion By Steel Fragments; Empiriaal Relationships for Fragment Residual 
Valoaity and Residual Weight," Technical Report #47, Ballistie Analysis 
Laboratory, Institute for Cooperative Research, The Johns Hopkins 
University, Aprtl 1961. 
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obliquity proposed by Otto P, Fuchs.? The remainder of this report is e 
a discussion of his model, the modification that has been made, and 4 
supportive arguments for use of the model. 4 
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II. THE RESISTIVE FORCE 


‘he expression for the resistive force acting during the penetration 
process as proposed by Fuchs involves a sum of three components that 
are functions of the instantaneous velocity: 


F=f, (v°) +f, wi) + £, (v2). (2) 


The first component, a static force, is defined as the product of 
the projectile cross-sectional] area and a stress factor. The stress 
factor is closely approximated by the target Brinell hardness when ex- 
pressed jn dynes/cm*. (The Brinell hardness number is multiplied by 
9.8 x 10’ to obtain the value in dynes/cm2). The first component is: 


ee ee eee 


TT ES 


0 ~ 
f, VW) =A. (2a) 


= 


The second component is a combination of the first and third 
components, hu.ce, the third component will be presented next. Analogous 
to the acrodynanis: resistive force, the third component is: 


. Bun 2 
f. (Vv) =CA Py vy ’ (2b) 
Sotuwtriint oo co the second component, it is defined to be: 
3 Bon O ta, uae, 

Poyty a % f 

five) #2 ye, ) £,(V") 
aa (2¢:) 
a atte pg See 
Ba = 7 1 & . 
oe zayc Ho, Vy 
Bee “his componenc can be compared to Stoke's Equation’ in which the resis- 

tive force is propertional to the velocity. 


A more general form for the resistive force equation has been 
adopted since the equation proposed by Fuchs did not yield satisfactory 
results when applied to available data. The generalized equaticn is: 


Mey A eT 


3Fuchs, Otto P., “Impact Phenomena," AIAA Journal, American Institute of 
Aeronauttes and Astronautics, Vol. 1, Nr. 9, Sept. 19638. 


“erence, M., Lemon, H., and Stephenson, R., "Analysical Experimental 
Physies," Universtty of Chicago Press, 1956. 
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PeA (C,H. +C, vit, 0, V+ Cy o, v; ). (3) 


The coefficients C,, Cz and C, will be determined in Section V. Equa- 
tion (3) can be expressed more succinctly as: 


_ etl 
Fe A (kK, + K, V+ K, V,"). (3a) 


III. THE DEPTH OF PENETRATION EQUATION 


Newton's second law of motion states: 


4 

_ dtm) _ av | 4 

p= Seo = uy +e ) 1 

‘ : ; 
assuming that the mass is held constant. The work done in penetrating i 
an increment "dx" is given by: 4 
are dV _ dx. ; 

Work = F dx = - mdx G7 =- mae dv -mv dv (5) 

‘ 

where the minus sign is due to deceleration, . 
4 

Substitution of Equation 3a for the force and arranging terms : 
yields: i 
> 4 

a vaVv at 

dx = =F TIRE eee RCT he (6) | F 

Ky + K, V+ K, \ <i 

4 

For a specified striking velocity and residual velocity, Fanation 6 is 
integrated to find the target plate thickness (ov the maximum depth otf ok 


penetration fur zero residual velocity). 


Vv 


x, Pe ‘ry 
Picea hs 
A : K, + Ky V+ K, Vv 


Ss 


Be 
de. 


ce einer Ba aS 


The psojectile mass and cross-sectional area are assumed to be constant 
in Equation 7. In most cases, there is very little, if any, mass loss 
when penetrating a single target plate. When deformation of the pro- 
jectile occurs, the cross-sectional area increases and should be 
accounted for. The effect of the cross-sectional area will be covered 
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in Section VII and will be assumed constant in performing the integra- 
tion of the equation. 


In order to integrate Equation 7, it is hecessgry to determine the 
value of the discriminant q, where q = 4 K, Kz - Kz. 


If q > 0, then: 


Xe Bee 5 a 3 f?Ks¥ + K\\I|“r 
X = ® (2) ln (K, + KV + KV’) - “Tz tan sa . (8) 
0 3 q q V 
s 
Substituting the limits and taking into account the negative sign 
yields: 
m K, + K.V_ + K,V 7 
acai 2 3's 
tc 2K, K, + KV. + KV. 
(8a) 
2K 2K,V_ + XK 2K,V. + K 
2 -1 3 r 2 -1 3s 2 
+ 1/2 tan ( 172 ). tan (Se) . 
q q q 
When q = 0 (which is the condition for Fuchs' original equation), 
integration of Equaticn 7 yields: 
Xe sae xi/2 V, 
X = ‘a (x) In («”? + re \ +5 1 3 ; (9) 
3 K + K V 
10 1 3 Vv 
s 
After substitution of the limits, Equation 9 becomes: 
Ki/2 & 2/2 V qi/2 xi/2 
1 $ 1 z 


Hee) PSII neater Ne Aone ten en 


, : 


Finally, if q < 0, Equation 7 is evaluated as: 


V 
2 (x) 2 we 2kV +k, - VAN} 
Xi. ae In OR + KY + V7) - (G2) in A Sen ees 
| A 2K. 1 2 3 /-a 2K,V + k, + 7G : 
0 Me 
(10) 
which is 
e 2 
: my 1 Ky + KV, + KV. 
Xe = (aK Po zZ 
; 3 Kee. RN. 
* (=) In Sie ae et Bele ak ae a oe) . (10a) 
=q (2K... + Ky + 7-q) (2K,V,, + Ky - 7) 
: IV. THE TIME-PENETRATION EQUATION 
The time to penetrate 2 target plate to a depth x can be found by 
: again considering Newton's second law of motion. 
oo ae (11) 


ae 


Substituting Equation 3a for the force, solving for the time, and 4 
: integrating (assuming that the projectile mass and cross-sectional area Bi 
; are constant) yields: 4 
: 
: : ee = ve q 
' i K, + K,V + K,V i 
0 Vv 1 2 3 
; a s 
f 


Bey ; 
hy { (the minus sign is due to deceleration). 
t 


Again letting q = 4k K, - coe three cases exist: 
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For q > 0: 


For q < 0: 
Tse 4 f ks"s +k, - ¥-q) (2K,V_ + K, + ¥=q) 
A \ rq LA Ce a) vq) (2K,V_* Ky = faq) 3 (15) 


whore T, is the time (in seconds) required for the fragment to penetrate 
to a depth x (or, in other words, until the velocity drops to the value 
V..). 

r 


V. DETERMINATION OF THE CONSTANTS C,, C. AND C 


1 2 3 

A non~linear least squares computer program was used to eveluate 
the constants C,, C, and Cz. The experimental data are used in 
Equations 8a, on or 10a (depending on the value of the discriminant, q) 
to detemnine the best values for C,, C, and Cx. To use the program, an 
initial guess is made for the values of the constants. The program 
computes new values for the constants based on values it computes for 
the partial derivatives of the plate thickness (the dependent variable) 
with respect to ths constants. The statements to evaluate the partial 
derivatives are provided to the program in 4 subroutine and are 
tabulated in this report in Table I{Ja. The computer program arrives 
at a convergent set of constants when the change in the value of each 
constant from one trial to the next becomes less than some predetermined 
tolerance value (0.01 was used in this case). 


Listed in Table IfIb is a summary of the computational runs made. 
Shown are the initial guess values for the constants with the corre- 
sponding root-mean square error, then the general set of constants 
obtained when combining the data for all target materials, and finally, 
the convergent set of constants obtained for cach target material. 

Also shown are the corresponding sigma and T-statistic test value for 
each constant of the convergent set. The data used are those tabulated 
in Appencix A with the exception of 20 out of the 277 detum sets. 
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Table IIIa. Statements To Evaluate Partial erivatives For 
Non-Linear Least Squares Program 


t Cc c 
C ¢ A = PROJECTILF CROSS=SECTIUINAL AREA (CiH#x2) 
¢ ¢, Cl = EMPIRICAL CONSTANT TO BE EVALUATED 
C s C2 s EMPIRICAL CONSTANT TO BE EVALUATED 
C ( C3 & EMPIRICAL CONSTANT TO BE EVALUATED 
or, SE i. HT & HARDNESS UF TARGET PLATE (OYNES/CMeeQ ) 
Co Me MP = MASS QF PROJECTILE (GRAMS) 
¢ ; Pl = FIRST OFRIVATIVE OF XT WITH RESPECT To C1 
re C P2 © FIRST VERTIVATIVE OF XT WITH RESPECT To C2 
Ee Pa = FIRST VERIVATIVE OF XT WITH RESPECT TO C3 
C ( RHOT « DENSITY NF TARGET PLalTE (GRAMS/CC) 
C ( VR « RESTHOUAL VELOCITY OF PROYECTILE (CM/SEC) 
C ( VS = STRIKING VELOCITY OF PROJECTILE (CM/SEC) 
C ' XT = TARGET PLATE THICKNESS (CM) 
C ¢ 
C C 
S2S)JKRT(RHOT#HT) 
yeSsayvys 
VERO T REVS em 2 
yeS*Vk 
LERHUTSRVRRE2 
1226045 
24 O8RHOTAHT 
42 Ska 
Re? O#eKHOTHYR 
T2MP/12.,08A*RHOT) 
S2.C*RHOTHYS 
Q0204,0%C14C 320 2H 42) eRHNT eHT 4 
Cc : 
1h (2C0.EG29.)GOTD 200 4 
¢ ¢ 
11 (wO.LT.LO.9GOTA 30N 4 
c - iw 
C ¢ m4 
C \, Thk DISCRIMINANT Q 75 GREATER THAN ZERO 
C f : 
VOC ILS) HAT OC 2HB4C aa i 
Yea PeHTo+C2e04+C aak te 
7325G6RT(Q0) i 4 
. J4sygeny e 
‘252° F¥R4C 205 } 
| Yost 2eeC Ques “4 
: Q724T0N(Q5/Q3) J 
{ NH2ATAN(Q6/Q3) 
| Q92T (03 
\ 19=C2/Q3 


QL1=C2#F#G6/Q4 
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Table IIIa. (Cont'd) Statements To Evaluate Partial Derivatives For 
Ron~-Linear Least Squares Program 


YL2-1.604(Q5/Q3) eae 
21321.94(Q6/Q3) e¥2 

W14205/Q4 . 

Q152=06/Q4 
T122%O114(Q8-07)/2,9 
T25Q)0Q1L1%(06/013°050129/2,.0 


pL eQO#(HT/OL“HT/C2*T AST?) 


Tax (F/Q34C 298 24F aH /Q4)%(Q798) 
143(5/034C28H¥014)/012 
T52(§/Q39C24H#Q15)/033 


223094 (B/QL=9D/Q247 34 QlLOwFS( 24-T5)) 


TASALTIG(Q2Z/Q1L)/C3 
T7=2C2"FH(C1HG6/ 07, 9%0G)4100/(9394C3) ) 21 Q5~07) 
T42(R/Q3—C146HQ14/2,0)/912 

192 (U/Q9=C146H15/2.0)/Q13 


PIV (C/AQLEVQZ4TOSTI+QLO¥F UCT R19) ) 


XT#I9¥(ALOG (21/02) 4C24F 199% ( 27-Q8)) 
RETURN 


THe DISCRIMINANT Q 1S EQUAL TO ZERO 


QwesSQKTCCL#HT) 
QLeQa+SQrT(CI#RHNTI*VS 
YZEQZ+SQRT(CI*RHNT) eR 
We VSeSQRTCHT/O1) 

Q5e: -SeSQORTCRHOT/CI) HVS 
Jee, 5S*SQkTLRHOT/CIISVR 
Qres Gat /C3 


PLEQTHIGH( 2g D/Q1 a2 eM /QZH+QI/Vc MH2Z—AF/QLHH2) 


Q4aceeS5QRT(HT/C3) 


925 248Q8/2 
Q1Q Q1L+«Q3 
Q11+Q2-Q3 
Qlz20"4+Q10 
Yi 3-Q94Q)1 


P22Q7HQ9H(2,0/Q12—2.0/Q13FQ98( 1 O/QILFF#241,0/Q91 20823) 
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Table IkIa. (Cont'd) Statements To Evaluate Partial Derivatives For 
, Non-Linear Least Squares Program 


TLRu3/0} "06/02 
12206 /Q24%9 695/91 482 
on € ; 
PISO7*(Q5/OL-G3/02"ALGG(G1/02)/C3+T14Q34T2) 
JF ak 
; XT 407 *(AL0G(Q1/02) 403/01 +Q3/02) 
RETURN 
{ THE GISCRIMINANT Q 1S LESS THAN ZERN 
f c 


ADO QTAIQRT(wQD) 
OLe LeHY SC 2eb4C aur 
Q22 LeHT oC aeO4C Jae 
327.086 SBR OC 245 a1)7 
My e342.0%Q7 
Q5e2, 080 SHU OC 285497 
Nes5=2,ONO7 
Ww 2S /Q7 
QIsT/C3 
QLOZALINGCQ3HQ5/(Q4*Q6)) 
QLL2C 36/02 ¢08Q7) 
Q12:C28H/07 
G1 32C146/62,0%Q7) 
FY eC BBOHOMQLO/ (~Q0)4Q19/034011/0Q4-Q11/95~0114/06 


C « 
PLedOeCHT/OL<HT/Q2+Q84T} ) 
¢ { 
Tees /Q lH CQHH/ (Q7KMS ) 
THe (S012) /Q94(54012) /Q8=(S4Q12) /Q4—(5~Q12)/Q6 
¢ ¢ 
PAPVORUBSQL§D/Q2+T2eQLOHQGHT A} 
a Gc 4 ' 
an Tr eQRBQLOR(LO/C RHC 1HG/ (2,0%ABS(Q0))) ‘ 
A Poe (ReQy 3) /Q34(UHQ19)/Q5 
he Tea (RmQ1 9) /044(U4Q13)/Q6 ; 
ras cof ) 
He! P32Q99 (C/QILME/Q2~ALNG(Q1/02) /C3=T4+Q8%( TS=TH) ) { 
a ' C ( \ 
lye XV a9e( ALONG (01/02) +Q8*Q1.) 4 
RETURN K 
i ENT ‘ 
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Table IIIb. Non-Linear Least Squares Fit to Thor Data 


Target Number of 
Material Initial General Convergent o _T _ Datum Sets 
Magnesium C, 0.40 0.70 1.96 1.72 1.1 22 
C, 0.90 0.23 “1.17 1.84 -0.6 
C, 0.80 0.50 0.83 0.42 2.0 ‘ 
ERMS 0.72 0.60 0.58 - - 
Aluminum C, 0.39 0.70 0.62 0.27 2.2 83 : 
C, 0.68 0.23 0.41 0.52 0.8 
C, 0.41 0.50 0.40 0.21 1.9 
ERMS 0.16 0.15 0.15 - . 
Titanium © C, 0.40 0.70 4.06 2.92 1.4 18 
C, 0.80 0.23 ~3.10 3.04 ~1,0 
C, 0.50 0.50 1.38 0.74 1.8 
ERMS 0.12 0.17 0.11 - - 
Cast Iron © C, 0.70 = 0.70 0.37 0.39 0.9 19 
C, 0.23 0.23 0.13 0.58 0.2 
Cz 0.50 0,50 0.53 0.18 0.3 
ERMS 0.18 0.18 0.09 - - 
eo ~ Steel (RHA) C, 0.40 0,70 0.53 0.93 06 17 
| C, 0.40 0.23 0.31 1.32 0.2 
Cz, 0.30 0,80 0.34 0.44 0.8 
ERMS 0.07 0,12 0.06 - - 
Steel (FHA) C, 0.00 0.70 -3.4x104 0.22 - .002 24 
C, 0.20 0,23 1.18 0.30 3.9 
7 C, 0.50 0.50 0.33 0.10 3.4 
ae ERMS 0.47 0.14 0.08 - - ; 
cial Copper C, 0.10 0.70 0.48 0.98 0.5 27 
me C, 1.50 0.23 0.20 0.68 0.3 " 
: Cz 0.50 0.50 0.61 0.09 6.7 
a ERMS 0.16 0.18 0.09 : - 
fie | 24 
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; Table IIIb. (Cont'd) Non-Linear Least Squares Fit to Thor Data 


Target Number of 
Material Initial General Convergent —o _T_ Datum Sets 
Lead Cc, -1.00 0.70 -7.70 1.84 -4.2 26 

C, 1.00 . 0.23 6.22 1.23 5.1 
Cc, 0.50 0.50 0.19 0.11 1.8 
ERMS 0.36 0.34 0.30 - - 
Tuballoy C, ~0.50 0.70 -0.34 0.59 -0.6 20 
C, 1.00 0.23 2.08 0.81 2.6 
C, 0.25 0.50 0.30 0.22 1.4 
ERMS 0.44 0.12 0.11 - - 
Combined — C, 0.40 0.70 0.70 0.12 537 257 
nese C, 0.90 0.23 0.23 0.13 1.8 
c. 0.50 0.50 0.50 0.03 18.6 
ERMS 0.29 0.23 0.23 - - 


Table IIIc. Summary of Convergent Values For 


C. C, and Cc, 
Target es °71 a3 "2 Sy °3 Nr. cS 
Magnesium 1.96 1.73 “1.17 1.84 0.83 0.42 22 0.58 
Aluminum 0.62 0.27 0.41 0.52 0.40 0.21 83 0.15 : 
Titanium 4.06 2.92 -3.10 3.04 1.38 0.74 #18 0.11 
Cast Iron 0.37 0.39 0.13 0.58 0.53 0.18 19 0.09 . 


Steel (RHA) 0.53 0.93 0.31 1.32 0.34 0.44 17 0.12 


Steel (FHA) -3.4x1074 0,22 1.18 0.30 0.33 0.10 24 0.08 
Copper 0.48 0.98 0.20 0.68 0.61 0.09 27 0.09 
Lead -7.70 1.84 6.22 1.23 0.19 0.11 26 0.30 
Tubal loy -0,34 0.59 2.08 0.81 0.30 0.22 20 0.13 


Combined Data 9.70 0.12 0.23 0.13 0.50 0.03 257 0.23 
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These 20 datum sets are tabulated in Table IV. 


Two things should be noted concerning the values of Table IIIb. 

First, the values for the constants can vary considerably with little 
change in the root-mean square error. Secondly, the sigmas for the con- 
stants are large relative to the value of the constants for the individ- 
ual target materials. This becomes more apparent in the summary Table 

‘ IIIc. Therefore, the constants evaluated with the combined set of data 
seem to be relatively good estimates for a general set of constants. 
These values are: C, = 0.70, Cz = 0.23 and Cz = 0.50. (The values for 
Fuchs' original equation are 1.0, 1.414, and 0.50 when the shape factor 
is unity. However, Fuchs' equation corresponds to Equation 9a rather 


than 8a). 


Substitution of the general set of constants into Equation 8a yields: 


2 
_ : m i 0.7 H, + 0.23 vii, P, V5 + 0.5 P, Vy 
t Ap 2 
t 0.7 H. + 0.23 Vi, P, V.. + 0.5 P, Vy 
— | 
ie pe, V. + 0.23 ) | 
£05306 | tan) Pe (16) | 
1.16 VH. p 
tt 
— ! 
_, [o, V. + 0.23 VAL p | 
ane ( pes t :] | 
1.16 vii, P. I. 
The penetration time equation becomes: [ 
| 3 
tN 
2m p, V_ + 0.23 VH, 9 ‘ :: 
0 il ae Puen Pigg ae tt 4 
. ‘i 1.16 “H ‘ 
1.16 A tt ‘ t Py 3 
an ; 


p, V_ + 0.23 VH. p 
- cant 7 :*s) : 


Equation 8a along with the values for the constants C,, C, and C, 
or in the form of Equation 16 will be referred to as the Z/F equation. 
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Table IV. Data Eliminated From Non-Linear Least Squares 
i Fit to Ther Data 
: Target Datum 
P Material Set Nr. He et M, Dd, Vs APs 
{(Kg/mm?) (cm) _{gms) (cm) m/s __ m/s 
; Magnesium 21 72.0 2.540 15.56 1.49 1417.0 537.7 
4 Aluminum 78 120.0 2.540 15.56 1.49 975.4 0.0 
b 
iM Titanium 1 190.0 9.127 1.95 0.76 $67.8 521.2 
" 1a 190.0 0.318 3.89 1.01 620.3 500.8 
Steel (FHA) 2 400.0 0.345 1,95 0.76 748.3 0.0 
4 409.0 0.345 1.95 0.76 1081.7 0.0 
9 400.0 0.635 3.89 1,01 1066.8 0.0 
14 400.0 1.270 3.89 1.01 1791.90 0.0 
Copper 18 42.0 0.318 7.78 1,27 745.5 501.7 
Lead 1 5.5 0.318 1.95 c.76 2401.2 1066.8 
2 5.5 0.318 1.95 0.76 2439.9 914.4 
4 5.5 0.348 1.95 0.76 957.7 457.2 
6 5.5 0.348 1,95 0.76 1721.5 762.0 
11 5.5 0.698 3.89 1.01 757.7 $33.4 
14 5.5 0.318 7.78 1,27 1810.5 1101.8 
25 5.5 0.635 15.56 1.49 2608.8 1005.8 
32 5.5 2.54 15.56 1.49 1263.7 762.0 
| Tuballoy 8 240.0 0.254 3.89 1.01 1471.6 823.0 
| 15 240.0 0.508 7.78 1.27 1699.9 1219.2 
22 240.0 0.508 30.16 1.74 2171.4 640.1 


Total number is 20 sets. 
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ae VI. COMPUTING RESIDUAL VELOCITY 


An iterative procedure must be adopted to sclve Equation 16 for re- 
Sidual velocity. One procedure is to determine the force acting ever con- 
secutive Ax incremerits and compute the corresponding speed reduction. 

Fron Equation 5, 


dv = ~ F dx/(m v) = - (Ky cn Kv") dx/(m v). 
The algorithm displayed in Figure 1 outlines the procedure. 


This model was chosen for computing the residual velocity because 
it allows flexibility in defining the force series There is no need 
to check the value of the discriminant q (q - KS } since iv does 
not appear explicitly in this approach. Aieiveclie. this model is an 
approximation to the integrated equation but little error is involved 
because of the imposed criterion of 


AV‘ - AV" 


Vu raviy | < 9-001. (18) 


Qne problem that has been encountered occurs when the residual 
velocity approaches zero. The required Ax increment to cause the quan- 
tity of Equation 18 tu ke less than 0.091 becomes progressively smaller. 
Hence, a lower limiting value must be imposed on Ax or on the residual 
velocity in order to terminate the loop cycle. A tabulation of tha 
computer program deck is given in Appendix C along with sampJe output. 


An alternative metho! is to use Equation 16 by progressively in- 
creasing or decreasing the iesiduai velocity valuc until the computed 
plate thickness yields the correct plate thickness to some degree of 
accuracy. A tabulation of such a program is included in Appendix C 
along with sample output of the program. 


Also in Appendix C are tabulccions of programs to find the plate 
thickness using Equation 82 and the Thor equation, and a program tor 
finding residual velocity using the Thor equation. Sample output is 
included for each program 


VIf£. VALIDATING ThE MODEL 


Qne method for demonstrating the accuracy of a proposed model is to 
plot the predicted residual velocity or the predicted plate thickness 
against the experimental value. A perfect prediction will lie on the 
diagonal line. Plots are presented in Appendix B comparing the Thor 
equation with the Z/f equation for residual velocity and for plate 
thickness. The Thor data tabulated in Appendix A is used for making 
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1. Compute the force Fy acting at depth x) based on velocity vy: 


2 
Fy = K, + KoVy + KeVy 


a 


2. Compute the decrease in velocity Av' based on Fis (F dx = mv dv): 


Av' = F_ Ax/m_v 
Pp 


1 1° 


3. Compute the force F2 acting at depth xz based cn the computed 


velocity at Xo» viz., Wwey- Avi: 


: 2 
Fy K + KV + Kav, ‘ 


4. Compute the decrease in velocity Av" based on Fo! 


4¢ 
Av"! = Fo4x/m, vas 


5. Test for |(dv' - av") /z (Av! + Av'"')| < 0,001. 


6. (a) Test condition satisfied, (b) Fails test. 
Set: V5 24; (av! + Av")/2, Set: Ax = Ax/2 
F = (Fy + Fi)/2, Returr to step 2. 
xX = xy + Ax, 


and repeat procedure unt: 1 
x = taiget plate thickness, 


Figure 1, Outline of Procedure te Find Residual Velocity 
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these comparisons. A study of the plots will show visually the accuracy 
and the similarity hetween the two modeis. 


The Thor equation plots do not necessarily look like a least sqvares 
fit to the data because the exponents were obtained from a least squares 
fit to all the data for each target material reported in Reference 2. 
These data includ? oblique angle turgets. What is plotted on the graphs 
of Appendix B js the data (tabulated in Appendix A) which involves 
normal impact only. 


The similarity in the predictions made by the two equations and the 
poor showing for some of the data for each target material is likely due 
to inaccuracies in the data. The striking velocity and the residual 
velecity can be in error because ct difficulties with the recording 
instrumentation. In some cases, the residual velocity was estimated 
from the depth of penetration into Celotex or similur material. Even in 
those cases where the residual velocity was determined from velocity 
screens and a chronograph, there can be doubt as to whether the same 
particle triggered both screens. 


A second parameter which is possibly inaccurate is the cross-~ 
sectional area of the projectile on impact. The yaw angle was not re- 
ported, resulting in uncertainty concerning the orientation of the pro- 
jectile. Table V lists the values for the cylindrical rod cross- 
sectional ereas as a function of yaw angle. For example, the cross- 
sectional area for the 1.95 gram cylinder at 0° yaw is 0.452 cm”, ‘he 
maximum area occurs at about 40° yaw and is 0.611 emf, at yo? yaw 
(sideways impact) the cross-sectional area is the mirimum ~ 0.411 cm*. 
The uncertainty in the area can represen: 24 wich es & 20 percent error 
since the value which was used in huti: the Thus equation and the 2/7 
equation is the value at 0° yaw. 


As indicated in Section 11, the cross-sectional area was assumed 
to be constant when pertorming the integration. In reality, the pro- 
jectile deforms and increases in cross-sectional area as it penetrates 
through the carget plate. However, in using the non-linear least squares 
program to evaluare the three constants, the effect of the increase in 
projectile cross-sectional area was statistically taken into account. 
In other words, the particular values of the constants which were 
selected represent an average effect of the projectile penetrating the 
target plate. 


The third parameter which is questionable is the target Brinell 
hardness. The values which were used are the nominal values reported in 
handbooks except an average value was used for cast iron and for the 
tuce-hardened steel Experience has shown that the actual Brinell 
hardness for a particular plate can vary bv at least 29% from the hand- 
book value and seems to be a function of plate thickness, at least in 
the case uf 2024T-3 aluminum and roiled homogeneous steel. 
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Table V. Cross-sectional Areas By Yaw Angel For Cylinder Rods 


Mass (gms) = 0.973 1.946 3.891 7.782 15.564 36.250 
Radius (c™s) = 0.296 0.380 0.506 9.633 0.746 0.372 
Length (cms) = 0,457 0.541 0.617 0.795 1.143 1.661 


Yaw Angle (degs.) Areas (cm?) 


0.275 0.452 0.804 1.259 1.748 2.389 
0.298 0.487 0.856 1.342 1.890 2.632 


10 0.318 0.517 0.901 1.414 2.018 2.856 
15 0.336 0.543 0.939 1.476 2.130 3.057 
20 0.351 0,566 0.969 1.527 2.226 3.236 
25. 0.364 0.584 0.993 1.566 2.305 3.389 
30 0.374 0.597 1.009 1.593 2.367 3.517 
35 0.381 0.606 1.017 1.608 2.410 3.618 
40 0.385 0.611 1.018 1.611 2.435 3.692 
4S 0.386 0.610 1.010 1.608 2.442 3.737 
50 0.384 0.605 0.995 1.580 2.430 3.755 
55 0.379 0.596 0.973 1,546 2.400 3.743 
60 0.372 0.582 0.943 1.801 2.351 3.703 
65 0.362 0.563 0.906 1.444 2.284 3.635 
70 0.348 0.541 0.862 1.376 2.200 3.539 
75 0.333 0.544 0.811 1.298 2.100 3.416 
80 0.314 0.483 0.755 1.210 1.983 3.267 
85 0.294 0.448 0.692 1.112 1.851 3.094 
90 0.271 0.411 0.624 1.006 1.705 2.897 
Maximum Error = +29% +26% +21% +22% +28% +36% 


- 1% -10% -29% -25% ~ 3% ~ 0% 


Area = nr? Cos a + 2rL Sin a 
where: r - radius (cms), 
L - length (cms), 
and a - yaw angle (degrees). 
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i Mathematical comparisons of the two equations are tabulated in 
Tables Vila, b, c and d. The definitions for the column headings are the 
following: 
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ae 
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Number: number (n) of datum sets for the target material, 


Mean: ‘the arithmetic average X = EX, /u, 


F Variance: (£(X;) - n(X)*)/(m - 1), 


Standard Deviation: = /Variance, 


SERRE BLS te SES REY 


D: the deviant (the difference between the predicted value and 
the experimental value), 


R. E.: the relative error (the deviant divided by the experimental 
value), 


and £ : denotes summation. 


Tables VIa and VIb are for plate thickness, and Tables VIc and VId are 
for residual velocity; the deviants and relative errors are presented 
respectively for both sets of tables. For some target materials, the 
Thor equation renders less error overall in its predicted values than the 
2/F equation. For other target materials, the Z/F equation is better 
than the Thor equation. In all cases, the two equations do not yield 
grossly different results from each other. 


The advantage of the Z/F equation is that it is more general in its 
application than the Thor equation. The Z/F equation may be used with 
some confidence for any case where the values of the parameters (target 
obliquity and Brinell hardness; projectile mass, cross-sectional area at 
impact, and striking velocity; and either plate thickness or projectile 
residual velocity) are known. By contrast, the Thor equation is limited 
to those projectile/target materials for which sufficient experimental 
data exists to evaluate the necessary empirical exponential constants. 
It is also limited to the range of values for each parameter for which 
data exists. As has been shown in Table II, the empirical constants 

vary from one target material to the next. The Thor equation is an 
expression involving parameters thought to be significant in the 
projectile-target interaction. On the other hand, the Z/F equation is 
based on an expression for the resistive force experienced by a pro- 
jectile while penetrating a target. The three constants which appear 

in the equation were determined by fitting the depth-of-penetration 
equation to experimental data using a non-linear least squares procedure. 
While these constants have been evaluated empirically, the same values 
are used for all the target materials. As a result, one feels more 
confident in applying the Z/F equation to projectile/target materials in 
general, 
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Table VIa. Comparison of The Mean, Variance, Standard 
Deviation, Sum and Sum of Squares of The Plate 
Thickness Deviants 


Target Standard 2 
Material ‘Number Mean Variance Deviation _£D =D 
Magnesium 23 0.129 0.753 0.868 2.966 16.952 Thor 
23 0.021 0.517 0.719 0.474 11.378 Z/F 
Aivminum 84 = -0.085 0.024 0.155 -7.127 2.606 Thor 
84 -0.066 0.016 0.129 = -5.533 1.737 Z/F 
Titanium 20 0.030 0.014 0.120 0.604 0.292 Thor 
20 0.033 0.024 0.155 0.665 0.480 Z/F 
Cast Iron 19 0.015 0.010 0.102 0.284 0.190 Thor 
19 -0.144 0.006 0.077 -2.738 0.501 2Z/F 
Steel (RHA) 17 -0.094 0.002 0.048 -1.598 0.187 Thor 
17 -0.085 0.005 0.070 -1.442 0.201 Z/F 
Steel (FHA) 29 0.012 0.015 0.124 0.363 0.436 Thor 
29 = =—-0.090 0.014 0.119 -2.598 0.630 2/F 
Copper 28 0.022 0.011 0.107 0.627 0.324 Thor 
28 0.118 0.013 0.116 3.314 0.754 2Z/F 
Lead 34 0,180 0.337 0.581 6.127 12.232 Thor 
34 -0.061 0.209 0.457 -2.085 7.034 2/F 
Tubal loy 23 0.010 0.049 0.222 0.237 1.084 Thor 
23 0.029 0.042 0.205 0.673 0.941 2Z/F 
Combined 277 0.009 0.124 0.352 2.483 34.303 Thor 
Data 277 = - 0.033 0.085 0.291 -9.270 23.656 2Z/F 
34 
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Table Vib. Comparison of The Mean, Variance, Standard 
Deviation, Swa and Sum of Squares of The Plate 
Thickness Relative Error 


Target Standard 2 
Material Number Mean Variance Deviation = RE. 2 (R.E.)” 
Magnesium 23 0.001 0.107 0.327 0.014 2.349 Thor 

23 0.003 0.115 0.339 0.061 2.531 2/F 
Aluminum 84 =-0.135 0.035 0.188 -11.323 4.470 Thor 
84 = - 0.113 0.026 0.160 -9,473 3.185 Z/F 
Titanium 20 0.024 0.072 0.268 0.473 -0.530 Thor 
20 = -0.027 0.086 0.293 -0.536 1.650 2Z/F 
Cast Iron . 19 0.020 0.013 0.115 0.371 0.244 Thor 
19 =-0.215 0.017 0.131 -4,090 1.189 2/F 
Steel (RHA) 17. --0.353 0.041 0.202 -5.996 2.770 Thor 
17. -0.238 0.023 0.150 -4.053 1.327. Z/F 


Steel (FHA) 29 0.024 0.048 0.220 0.684 1.366 Thor 
29 =~ 0.083 0.019 0.138 -2.416 0.738 2/F 


Copper 28 0.017 0.033 0.182 0.483 0.900 Thor 
28 0.208 0.029 0.170 5.836 2.000 2Z/F 
Lead 34 0.307 0.599 0.774 10.428 22.964 Thor 
34 0.100 0.218 0.467 3.393 7.549 2/F 
Tuballoy 23 0.034 0.217 0.466 0.790 4.798 Thor 
23 0.074 0.186 0.431 1.693 4.212 z/F 
Combined 277.2 -0.015 0.149 0.386 -4,076 41.232 Thor 
Data 277 -0.034 0.087 0.295 ~9.579 24.381 z/F 
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Table Vic. Comparison of The Mean, Variance, Stundard 
Deviation, Sum and Sum of Squares of The 
Residual Velocity Deviants 

Target Standard 2 

Material Number Mean Variance Deviation ED tDY 

Magnesium 23 34.4 73596.7 271.3 790.8 1646313.1 Thor 


23 14.6 38264,1 195.6 336.9 846745.9 2Z/F 


Aluminun 84 -51.9 4577.3 67.7 -4361.4 606365.9 Thor 
84 -39.1 5763.3 75.9 ~3286.9 606962.7 2Z/F 


Titanium 20 42.8  20346.8 142.6 855.5 423167.2 Thor 
20 33.3 23154.6 152.2 665.7 462092.6 2Z/F 


Cast Iron . 19 -0.4 8712.5 93.3 -7.3 156828.0 Thor 
19 =-149.8 11266.8 106.1 -2846.8 629338,2 2Z/F 


Steel (RHA) 17s -119.9 4466.9 66.8 -2039.1 316062.1 Thor 
17, -100.1 5978.1 77.3  =-1701.5 265952.4 Z/F 


Steel (FHA) 29 -0.7  17081.4 130.7 “19.5 478292.9 Thor 
29 ~94.3  37147.3 192.7 -2734.9 1298040.9 2/F 
Copper 28 18.2 12024.6 109.7 510.7 333979.7 Thor 
28 89.2 14203.1 119.2 2497.2 606139.8 2/F 
Lead 34 82.7 60913.1 246.8 2813.4  2224937.2 Thor 
34 34.8  66762.6 258.4 1184.0 2244396.9 2/F 
Tubal loy 23 -23.6 59978.1 244.9 -543.3 1332350.9 Thor 
23 12.7) 72507.2 269.3 291.8  1598861.4 2/F 
Combined 277 -7,.2 27187.9 164.9 -2000.4 7518297.0 Thor 
Data 277 -20.2 30600.0 174.9 -5594.5 8558583.8 2/F 
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Table VId. Comparison of The Mean, Variance, Standard 
Deviation, Swn and Sum of Squares of The 
Residual Velocity Relative Error 


Target Standard 2 
Material Number Mean Variance Deviation ZR, BE. & (R. E. 
’ Magnesium 23 0.048 0.383 0.619 1.111 8.482 Thor 


23 0.062 0.259 0.508 1.431 5.778 2/F 


Aluminum 83* -0.113 0.043 0.207 ~9.370 4.577 Thor 
83* -0.085 0.073 0.270 -7.062 6.596 2/F 


Titanium 20 0.063 0.066 0.258 1.261 1.341 Thor 
20 0.059 0.103 0.321 1.173 2.032 z/F 


‘ Cast Iron 19 0.003 0.044 0.210 ~9,.049 0.795 Thor 
19 -0.352 0.041 0.376 -6,.638 4.901 Z/F 


Steel (RHA) 17 -0.189 0.015 0.121 ~3.214 0,844 Thor 
17. -0.139 0,011 0.103 ~2.366 0.500 2/¥ 


Steel (FHA) 26* = -0.049 


0. 0.149 ~1,.286 0.615 Thor 
28* -0.237 0. 
0 


22 
39 0.373 -6,.641 5.239 Z/F 


Copper 28 0.025 6021 0.144 0.696 0.576 Thor 
28 0.171 0.031 0.175 4.789 1.646 2/F 


Lead 34 0.166 0.150 0.387 5.664 5.880 Thor 
34 -0.001 0.126 0.355 -0,023 4,156 2/F 


Tuballoy 23 0.023 0.100 0.316 0.530 2.215 Thor 
23 0.025 0.173 0.416 0.578 3.286 2/F 


Conibined 273* ~-0.017 0.093 0.305 -4,657 = 25.325 Thor 
Data 275* 0.054 0.124 0.352 -14.809 34.764 2/F 


w 
Does not include data where the residual velocity was zero unless the 
predicted residual velocity was also zero. 
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VIIL£. ANALYZING THE PENETRATION PROCESS 


Altkough the force Equation 3 may not represent reality, it does 
| involve relevant physical parameters. When applied to Newton's equation 
of mction, the resulting integrated equation doss a fair job of predict- 
| ing piate thickness (1T residual velocity as the case may be) over a 
wide range of target materials, It would be fair to say that the force 
equation represents an average effect of the resistance to penetration 
encountered by a projectile. Therefore, it should be possible to learn 
something about the penetration process. 


The equations of Section III can be reduced to the following form: 


XA 


t 
ca = £(p,, He» Vo» V,) - 


t?’ 

5 The parameters associated with the projectile (the cross-sectional area, 
A, and the ‘projectile mass M,) appear on one side of the equation along 
with the target plate thickness X,. This allows generalized curves to 
be drawn, plotting the values for tee as a function of velocity for 
particular values of target plate hardness H, and target density Pps 


A family of curves is shown in Figure 2a. Tne target plate hardness 
is held constant at BHN 100. Each curve represents a different target 
plate density. Along the abscissa (the x-axis) is the parameter Ke A/Mp 
and along the ordinant (the y-axis) is the projectile velocity. The 
plot can be used in two ways. 


First, an estimate of the residual velocity for a given striking 
velocity can be made by finding the point on the curve corresponding to 
the striking velocity, then following the curve for a distance in the x 
direction corresponding to the computed value for X,A/M, and then read- 
ing the residual velocity off of the y-axis. For example: 


Given: xX. “ 1.5 cn. 

H, = BHN 100 = 9.8 x 10” dynes/cm 

A= 2.5 cm? 

Vs = 2500 m/s = 250000 cm/sec 

e.* 2.77 g/cc (aluminum) 

M = 2.2 grams. 

5) (2.5 
Find the residual velocity where X,A/Mp « AL 5 - = 1,70 by 


referring to Figure 3a. 


Proceeding over 1.70 units in the x-direction along the curve for 
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aluminum (BHN 100) beginning at the point for 2500 m/s yields a value 

for velocity which is less than zero. That is, the projectiie would not 
be able to penetrate 1.5 cm. Changing X; to 0.5 cm yields a new value of 
0.568. Using this value results in an approximate vusiduai velocity of 
Vy = 780 m/s. 


The second way these curves can be utilized is to estamate the limit 
velocity, i.e., the striking velocity required for the projectile to 
travel completely through the target with zero residual velocity. For : 
the second case in the example above, where X,A/M, = 0.568, the limit 
velocity is found by going to the left 0.568 units along the x-axis from 
where the curve for aluminum meets the axis and finding the corresponding 
x-point on the curve and then reading the velocity value on the y~axis. 
In this example, the limit velocity is approximately 1630 m/s. Greater 
precision (mathematically) can be obtained by substituting the appro- 
priate values into Equation 16. 


It should be noted in Figure 2a, Figure 2b (which is similar to 2a 
except the hardness is held constant at BHN 200) and in Figures 3a, b 
and c that the slope of each curve becomes more negative as the velocity 
decreases and changes dramatically as the velocity approaches zero. The 
latter is caused by the static component of the force equation dominat- 
ing the other two terms. This explains why it is not possible to 
linearly extrapolate limit velocities from residual velocity data. 


It should also be noted that the variation of penetration with 
respect to target plate density is non-linear. There is greater varia- 
tion fox a change in less dense materials than a corresponding change 
in the more dense materials. For example, a velocity change from 2500 
m/s to 2000 m/s yields the following: 


Density XpA/M, Pe X-A/M, 
1 g/ec 0.300 0.300 
2 g/cc 0.175 0.350 
4 g/cc 0.095 0.380 
8 g/cc 0.052 0.416. 


If the variation were linear, the values in the third column would be 
identical. A non-linear effect can also be seen in Figures 3a, b and c 
with respect to target plate hardness for a given plate density. Hence, 
the homogeneity of the target with respect to hardness and density 
determines to a great extent the replicability of a given set of 
experimental conditions. 


The single most important geometric variable of the projectile 
parameters is the projected cross-sectional area of the projectile at 
impact. The yaw angle is difficult to control and to determine with 
the excepticn of spheres and spin stabilized projectiles. Any change 
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in the yaw angle results in a new cross-sectional area since it is the 
projected area of the projectile onto the surface of the target plate 
that is required. 


APTS ere 
7 


When consideriig the ability of the projectile, as a whole, to 
penetrate a target, the important parameter is the ratio of the pro- 
jectile mass to its area. Shaped charges are capable of deep penetra- 
tions because they have a large mass per unit area ratio. 
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IX. SUMMARY AND FUTURE AREAS OF INVESTIGATION 


An analytic model of kinetic energy round penetration has been 1 
presented. This model compares favorably with the Thor equation in { 
predicting residual velocity for a projectile-target interaction. It is : 
more general in its application than the Thor equation and can be used " 
to study the penetration process. i 


An extension will be made to include oblique attack angles. A 
preliminary approach will be to adopt the same method as the Thor 
approach, i.e,, multiplying by the secant of the angle raised to some 
power (sec 6). Something more compiex may be required to adequately 
predict the effect of oblique angles. 


“AIS ne 


Te 


Another area of investigation is predicting projectile breakup and 
predicting residual masses. This information can then be used to pre- 
dict penetration of a secondary target plate. 
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APPENDIX A 


PROJECTILE-TARGET PENETRATION DATA 


All values reported in Reference 2 appearing in this appendix have 
been converted from the British system of units to the metric system 
with the exception of the Brinell hardness numbers which were already in 

: the metric system. 
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it te cme merc EE EME YA Om RS Ce ML AMAL 6 su T Amie NaN aa 


TH ARAR AO AOAOADKAARANKHAANANGAAAANKHAAAANAGAMAANANAARAGHAAGANAANH 


1 

2 

: : 3 

2/F FORCE PENETRATION MODEL : 

6 

GLOSSARY l 

PTT] IOENTIFIES REQUIRED INPUT DATA 9 

* IDENTIFIES INPUT DATA WHICH IS NOT REQUIRED He 

12 

13 

A = Frrst TERM GF THE FORCE EQUATION (K1) 14 

AA = STATIC FORCE COMPONENT (DYNES) 15 

eee ANG = IMPACT ANGLE (DEGREES) J6 

AME - CONSTANT BASED ON LEAST SQUARE FIT TO DATA 17 

B = COEFFICIENT OF SECOND TERM UF THE FORCE EQUATION (K2) 18 

aMF = CONSTANT BASED ON LEAST SQUARE FIT TO DATA 19 

‘ = COEFFICIENT OF THIRD TERM OF THE FORCE EQUATION (K3) 20 

CME ~- CONSTANT BASED ON LEAST. SQUARE FIT TO DATA 21 

vHELT ~ DELTA TIMF INCREMENT 22 

CELTAV = DELTA sPEED DECREMENT 23 

DELY = COMPUTED DROP IN SPEED BASED ON VR1 AND Fy 24 

GELVP - COMPUTED DROP IN SPEED BASED ON VR2 AND F2 25 

CELX - X INCREMENTAL VALUE 26 

DEVANT = THE COMPUTED RESIDUAL VELOCITY MINUS THE EXPERIMENTAL 27 

RESIDUAL VELOCITY ( THE DEVIANT ) 28 

+) = FORCE ACTING AT START QF X INCREMENT 29 

62 = FORCE ACTING AT END OF X INCREMENT 30 

VORCE = TOTAL FORCE 31 

VI - COMPONENT OF FORCE PROPORTIONAL Tu SPEED 32 

V2 = COMPONENT OF FORCE PROPIIRTIUNAL TO SPEED SQUARED 33 

* DTA - HOLE DYAMETER (CM) 34 

ICNT ~ INDEX COUNTER 3b 

* IR ~ SHOT IDEN IFICATION NUMBER 36 

PCA ~ PROJECTILE CROSS@SECTIONAL AREA (SQ CM) 37 

eee OP OTA ~ PROJECTILE DIAMETER (CM) 38 

kee LSS ~ PROJECTILE MASS (GRAMS) 39 

PRESUR = STATIC PRESSURE ACTING ON PROJECTILE (MEGA»PASCALS) 40 

RATIO - THE COMPUTED RELATIVE ERROR 4\ 

& RMASS = PRIMARY RESIDUAL FRAGMENT MASS (GRAMS) 42 

SUMET - TOTAL ¢MPULSE TO DEPTH X (UYNE~SECS) 43 

wok T BRN ~ TARGET PLATE SRINELL HARONESS NUMBER (KG/MM*¥#2) 44 

rEST ~ VALUE TO BE TESTED 45 

ee THICK  YARGET PLATE THICKNESS (CM) 46 
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Tse mene - 


AAMAAANQANAAARAANGA 


ana 


*he 


xe 


af 


cle 


9S 


L1O0Q 


110 


— 


PENETRATION TIME (MICRU SECONDS) 


TIME - 

TRHO - TARGET PLATE DENSITY (G/CC) 

VR = RESIDUAL SPEED (CM/SEC) 

VRY ~"RESIOUAL SPEED AT START OF X INCREMENT 
VR2 ~- RESIDUAL SPEED AT END OF X INCREMENT 
VRE = EXPERIMENTAL RESIDUAL SPERD (M/S) 
VRP ~ RESIDUAL SPEED (M/S) 

VS ~ STRIKING SPEED (CM/S) 

VSP * STRIKING SPEED (M/S) 

x = DEPTH OF PENETRATION (CM) 

XPRT ~ DETERMINES WHEN TO PRINT QUTPUT 
XPRT] - 


THE PRINT INCREMENTAL VALUE 


DATA P1/30141592654/ 
AMF 20,70 

BMF =0,23 

CMF 20,50 

WRITE( 6260) 


FURMATCLHL/LHOs20X240HSAMPLE OUTPUT FUR Z2/F FORCE PENETRATION » 


35H MODEL PREDICTING RESINUAL VELOCITY /) 
CONTINUE 
READ CARD WHICH ITDENTIFIES TARGET MATERIAL 
REAL (CS, 9SYTGTLaTGT2s TGT3 
WRITE(Os95)TGTIsTGT2sTGT3 
EMRiMAT(3A6) 
TONTSO 
CONTINUE 
A BLANK CARN SEPARATES TARGET MATERIAL GROUPS 
ENN PUNGHED IN COLUMNS 1 TO 3 WILL TERMINATE THE PROGRAM 


REOO (bs LLOYNR» TRHOs TRHNG THICK» ANGs PMASSs PUTA, VSPsVREsRMASSsHIIA 


EIRMAT( ABS FOe2aF Og lak, BaF Dela F Oy laF Oe de ZF, le 2F7,2) 
tF(TBHNeLE.O.NIGOTU 900 
GUNVERT STRIKING SPEER TO CM/SEC 
VS=VSP*100.0 4 
SOMPUTE CRUSS@SECTIONAL AREA 
PCAPL#(PDIA/2¢0) 4*2 
COMPUTE Kl» K2 AND K3. ¢ # AsBal) 
429.8 74TBHN 
A=SORT(CA®TRHO) *BMF 
Cs TRHUSCME 
as A*ANF 
AAsPCARA 
VRevS 
DETERMINE PRINT FREQUENCY AND INITIAL DELTA X INCREMENT 
XPRIT=THICK/1000 
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ae 


aoa 


DELAsTHICK/10,0 
IF(DELX .GT. 0201) DELX#Q.01 
INITIALIZE VARTABLES 
Xs0.0 : 
XPRT#O,0 
TIME*s0.0 
SUMFT#0,0 
TCNT@ I CNT 41 
WRITECOsS LIS) ICNT oNR 
115 pReNy oe oc AngescnesanORtUs SET NR» 25a 10Xs L2HIDENTIFJER ©92X%2A6/) 
WRITE (60180) 
150 FORMAT (1HOs9Xs LHX» 8XaZHVRe LOXSSHF (VO) s LOXe SHE (VL) » LOXe SHE (V2) 9 10X— 
} SHFORCEs OX, 4HTIMEs 4X, B8HPRESSURESSXs7HIMPULSE ) 
WRITE (69151) 
151 FORMAT(CLM 96Xs4H(CM) 2 5X2 5SH(M/S) 2 OXs THC DYNES ) » 8X0 THC DYNES) 2 8X. 
1 THC DYNES) 9 8X» 7TH( DYNES) 2 4Xs BM(MUMSEC) 26Xe18H(M PA) (DYNE~SEC)/) 
COMPUTE FORCE AND STATIC PRESSURE ACTING AT X20 BASED ON VR*VS 
FVLaPCASB¥VS 
FV2=PCAMCRYVS§ a2 
FORCE =AASFVL +FV2 
PRESUR@FORCE/PCAS) Eo? 
WRITE (69251) X,VSPyAA,FVIZFV25 FORCE, TIMEsPRESURs SUMFT 
200 CONTINUE 
IF(X .GE. THICK) GOTO 300 
TFCVR.LT.39e0)G60TO 290 
THE POLLOWING STATEMENT INSURES THAT THE PLATE THICKNESS 
WILL NOY BE EXCEEDED 
IFC CX#0ELX) .GT, THICK) DELXsTHICKeX 
IF(UELX «LEs 060) GNTO 300 
COMPUTE FORCE ACTING AT Xi BASED ON V1 
VRLEVR 
FL*PCA®( AFBMVRI+C UVR) OH?) . 
270 CUNTINUE 
COMPUTE DELTA v' BASED ON Fi 
NFLVSELEDELX/(PMASS#YR)) 
VR23VR=DELV 
TF(VR26472020)G60TO 290 
COMPUTE FORCE F2 ACTING AT X2 BASED ON V2"VL-DELV! 
F2aPCAM(BTAMVRZ4CHVR2 REQ) 
“OMPUTE DELTA v'!' BASED ON F2 
VELVP=F2eDELX/(PMASS&VR2) 
MAKE TEST 
TESTE(DELV=DELVP)/ ((DELV+DELVP)/2,0) 
LF(\BS( TEST) .47.20.001)G0TN 250 
TE(LEULX LT LeE*5)GOTO 280 
FAILS TEST © REOUCE DELX BY HALF (MINIMUM 1§ 0.00001 CM, ) 
NELXDELX/2,0 
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C 
C 
C 


250 
281 


260 


240 


aC 


GUTO 220 


CON) INUE 
FORMAT(LM 9F10e3sF10,1,1P4E15-520PF12.3,2F12.1) 
MEETS TEST - 


DELTAV#(DELV*#DELVP)/2e0 
COMPUTE VReVI-DELTA V 
VRaVR-DELTAV 
FIND THE AVERAGE FORCE ACTING OVER THE SNTERVAL 
FORCE=(F1%F2)/2.0 
UPDATE THE DEPTH GF PENETRATION X 
XaX+DELX 
UPCATE THE PENETRATION TIME: IN MICRO-SECONDS 
NOEL T=PMASSH#DELTAV/FORCE#].E6 
TIME @TIME*DELT 
UPDATE THE ToTAL MPULSE EXPERIENCED BY THE PROJECTILE 
SUMFT®SUMFY*FORCESDELT #1.E-6 
CHECK FOR PRINTING INFORMATION 
IF CX ,GE.XPRT)GOTO 260 
GUTN 200 
CONTINUE 
RESOLVE THe VARIARLE COMPONENTS OF THE FORCE 
FVLaPCA*BY(VRODELTAV/2.0) 
FV2=PCASCM(VR4+DELTAV/2,0) ¥¥2 
OMPUTE THE STATIC PRESSURE 
PRESUR#BFORCE/PCAW] Ew? 
VRP=EVR/100.0 
WRITE (62251) Xa VRP SAAS FVL es FV2,>FORCEs TIME, PRESUR»s SUMFT 
XPRTSEXPRT+XPRTI 
SOMTINUE CYCLING UNTIL TARGET PLATE IS COMPLETELY PENETRATED 
OR UNTIL THE RESIDUAL SPEED IS LESS THAN 10 CM/SECe 
GOTU <O0 
CONT IAUE 
VETAULT = SINCE THE RESITOUAL SPEED IS LESS THAN 10 CM/SEC 
VReEO.U 
DELTAVe 1.9 
CONTINGE 
COSPUTE AND PRINT FINAL VALUES 
FVL=PCARBR(VREDELTAV/2,0) 
FV2=PCAMC#(VRODELTAV/2,0) #2 
SUMF TeSUMFT*+FORCERVELT 1.64 
PRESURZFORCE/PCAM)] eEw? 
VRP -VR/100.0 
WRIIE(Os251) Xs VRPZAASFYV] »FV2, FORCE, TIME» PRESUR»s SUMFT 
CONTINUE 
UMPUTE AND PRINT THE DIFFERENCE BETWEEN THE PREDICTED VALUE 
“NU THE EXPERIMENTAL VALUE ( DEVIANT ERROR) AND THE RELATIVE 
tRROR 
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14) 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
\77 
178 
179 
180 
181 
162 
183 
184 
185 
186 
187 


+s 


- 


nanan renal Smut ey 


DEVANTaVAP@ VRE 
RAT 10"1030.0 
IFCVRE LE. 1.9) GOTO 302 
RAT LOMDEVANT/YRE 
302 CONTINUE 
WAI TE (65310) THICKs VRESDEVANTS RATIO 
BLO FORMAT(LHOsPLOs Se FAC, L23Xe 22H" EXPERIMENTAL VALUES » 8X» 
1 17HSPEEO DEVIANT 1s pELOcL*2Xe23HRELATIVE SPEED ERROR 18-F10,3) 
PRINT THE INPUT VALUES 
WRITE (6,328) 
215 FORMAT  LHOs LOXs HT SHNs SXs SHTHICKs 3X» SHPMASS » 7Xs BMANGs TX» SHVSP 
} 7TXs SHVRE» SXo 5SHRMA 5S, 6X2 4HPDI As OX» 4HHD IAs OXs 4HTRHO ) 
WRYTE (09920) TRAN, THICK, PMASS» ANG» VSP» VRE, RMASSJ POL ANDI Ay TRHD 
320 FORMATO()H slOHINPUT 1S = sFlOeds2PlOedeFLO Ls 2FLO.bedF10,32F10,2) 
CONTINUE WITH NEXT CASE 
LFCHODCIOCNTs 2) eNE.OIGOTO 100 
WRITE (6960) 
WRITE (Gs 98)TUTIsTGT2, TTS 
GNTO 100 
R00 IF(NR wNEs OHEND 1GO070 90 
WRITE (69901) 
90) BURMAT(LHOs20Xe LOHEND OF RUN ) 


STOP 
END 
DATA 
TANTUM ee 
160 4.46 0.0 ie? 20 «1,95 6759 567,864 $821.21 1.88 
4.0 4,48 190.0 2348 Oo 1,93 “$9 880,26 390.40 ©10,09 
he 4.44 190.0 oak 10 #1695) ./59 1335,65 1063.56 1.63 


0 
He 4,48 190,90 2695 -O 1.95 «739 1491,08 683.36 1073 

S 4,48 140.0 2638 09 =1,95 6759 1986,38 1127.18 056 
7.0 4 4b 190.0 2 2371,65 $81.61 10.00 
9 4,48 190.0 0127 -O 3.89 1,013 798,866 672,69 3.85 
9 4.44 190 9 «127 -9 3,89 1,013 1032.66 937,68 Lee 
O 4.48 190.0 2318 09 = 3,89 1.018 £20.27 300,79 3,43 
1 4.48 190.0 2318 . 3.8% 1.013 773,20 $92.78 3,81 
q 
ia) 
9) 
\) 


. 
nn 
oS 
ea 
o 
~ 
‘7 
bd 
wu 
eo 
=~ 
uw 
<s 


9) 
4.48 190,9 2 aha O 3,89 16013 1499,01 1251.81 2043 
4.44 190.0 2635 © 3,89 16013 1505,71 774.50 3424 
4048 190.0 +35 9 3.89 16013 1526,13 979.32 ~10.00 
4.4A 190.0 0695 9 3,89 16013 2455.16 1367.33 oe 
16.9 4.48 190.0 1.270 0 3,89 1,013 2551,79 1165.86 057 
i7e 448 190.0 127 O 7.78 16267 640,99 361.44 Tez 
O 7.78 16267 959,21 874,78 7.72 
O 7,78 16267 975,97 785.47 7672 
Q 7.78 12267 1484,38 994.87 ~10.00 


LHe =4.48 190.0 » 127 
160 4.48 190.0 2318 
21049 4.48 190.6 695 


or 2s > @ = @ ww @ 


BY 


RHA 


7.78 
7.78 
7.78 
7.78 
7.78 
7.78 
7.78 
7.78 
7.78 
7,78 
7.78 
7.78 
7.78 
7.78 
7.78 
7.78 
7.78 


135.0 
133.0 
300.0 
305.6 
393.0 
133.0 
133.0 
135.0 
300.0 
300.0 
300.0 
300.0 
300.0 
300.0 
300,0 
305.0 
332.0 


1.95 
1.95 


1,93. 


1.95 
1.93 
3.69 
3,89 
3,39 
3.89 
3,89 
3.89 
7.78 
15,56 
15.56 
15,56 
15,56 
15.56 


90 


0759 

e759 

2759 

0759 

0759 
1.033 
12013 
12013 
16012 
16013 
1.013 
1.267 
1249) 
1249) 
1249) 
1249} 
1249} 


688.49 
1211.58 
1521,26 
1394,46 

609,90 

302,06 

393,50 

883,62 

879,65 
1466,09 
1466.09 

909,52 

916,53 
1425,55 
1432.56 
1556.00 
1660,55 


846.87 
1014.98 
1196.34 

460.25 

367.89 

277.37 

256.03 

$42.54 

563.69 
1166.34 

687.32 

690.37 

754.38 
1179.58 
1037.8% 
1109.47 

759.56 


=10.09 
~10,00 
#10.00 
710,00 
~10.00 
#10,00 
#10,00 
“10.00 
#10.00 
#10.90 
#10.00 
#10,00 
#10.00 
©10,00 
10.00 
710.00 
#10.00 
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PROGRAM USING 2/6 EQUATION TO PREDICT RESIOUAL VELOCITY 


OLOSSARY OF VARIABLES 


IDENTIFIES REQUIRED INPUT DATA 
IDENTIFIES INPUT DATA WHICH JS NOT REQUIRED 


VBR 
uEL TAX 
ELV 
$b 
THR 
NVAR 
EBAR 
egh 
TVAR 
HOTA 
tcT 
oma) 
tClE 
{DW 
1FLGM 
TFL GP 
TREP 
POIA 
PMASS 
RATIN 


RELERR 
RMESS 
SUMS 
Surnse 
SUPRE 


praeouestoeerseeeratae ete tava eee 


THE ANGLE OF THE TARGET PLATE WITH RESPECT TO 
LINE OF FLIGHT = OBLIQUITY (DEGREES) . 
THE PRUJECTED CROSS=SECTSUNAL AREA OF PRGOVECTILE 
QN IMPACT : 
CONSTANT BASED ON LEAST SQUARE FIT fO THOR DATA 
SONSTANT BASED ON LEAST SQUARE FIT TO THOR GATA 
CONSTANT BASED ON LEAST SQUARE FIT TO THOR DATA 
THE NEVIANT = COMPUTED VALUE MINUS EXPERIMENTAL VALUE 
THE AVERAGE (MEAN) VALUE UF THE DEVIANTS 
ALLOWED TOLERANCE ON THE PLATE THICKNESS (CM) 
THE INCREMENT ON THE RESIDUAL VELOCITY 
THE STANDARD DEVIATION OF THE DEVIANTS 
CONVERSION FACTOR + DEGREES TO RADIANS 
THE VARIANCE UF THE DEVIANTS 
THE AVERAGE (MEAN) VALUE OF THE RELATIVE ERROR 
THE STANDARK DEVIATION OF THE RELATIVE ERROR 
THE VARTANCE OF THE RELATIVE ERROR 
THE oLAMETER OF THE HOLE MADE IN THE TARGET (C*) 
INDEX TO COUNT NUMBER OF DATA CARDS FOR ONE TARGET 
}NOEX COUNTER ON NUMBER OF PCIINTS FOR OEVUANTS 
INDEX COUNTFR ON NUMBER OF PLINTS FOR RELATIVE FRRUR 
AN IDENTIFICATION NUMBER GR SYMBOL=DESIGNATES SHOT NRe 
FLAG TO INDICATE RELATIONSHIP ON THICKNESS FOR THIS VR 
FLAG TO {NDICATE RELATIONSHIP ON THICKNESS FOR THIS VR 
NUMBER OF ITERATIONS COUNIER INDEX 
DIAMETER OF THE PROJECTILE (CM) 
MASS OF PROJECTILE (GRAMS) 
USED TO DETERMINE A FIRST ESTIMATE OF RESIDUAL 
VELOCITY 
THE RELATIVE ERROR OF COMPUTED VS. EXPERIMENTAL 
THE RECOVERED PROJECTILE MASS (GRAMS) 
THE SUM OF THE DEVIANTS 
THE SUM MF THE DEVIANT SQUARED 
THE §UM OF THE RELATIVE SRROR 
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SUMRES = THE SUM OF THE RELATIVE §RROR SQUARED 
wee = «6 TBHN)— om THE TARGET BRINELL HARDNESS NUMBER (KG/SQ@ MM) 
wee = =6THICK = THE THICKNESS OF THE TARGET PLATE (CM) 
whe TRHO = THE pENSTTY OF THE TARGET PLATE (G/¢C) 
VR - PREDICTED RESIDUAL VELOCETY (M/S) 
* VRE = THE EXPERIMENTAL RESJDUAL VELOCITY (M/S) 
aay VSP = THE EXPERIMENTAL STRIKING VELOCITY (M/S) 
xT ~ PREDICTED VALUE OF PLATE THICKNESS FOR CURRENT VALUE 


OF RESIDUAL VELOCITY 


XTO ~ VALUE OF PLATE THICKNESS WHEN VR®O (MAX PENETRATION) 


DATA PLsDOTOR / 3.1415926542 060174534925 / 
10 FORMAT( Abs2FoelsFO.3rF5sle2FOedrAFOela2F7,2) 
12 FURMAT(10A6) 
C1#0.70 
6230.23 
C3*s9,50 
DELTAX20,00] 
100 CONTINUE 
WRITE (69103) 
103 FORMAT(LHL) 
WRITE(62106) 
106 FORMAT(1HOs20X242HSAMPLE A FOR 2/F EQUATION PREDICTING » 
1 16H RESIDUAL VELOCITY / 
READ CARD DESCRIBING THE TARGET MATERJAL 
REA! (Ss L20TG6T1LaTGT22TGT3 
WRITE( Gs 12)TGT1lsTGT2,TGT3 
WRITE( 69110) 
110 FORMAT(1HOs4X%s39HNR TRHO TBHN THICK ANG MASS PDIA »6X> 
1 ZHVSa2X%eL2HRMASS HOTA 9 5X2 3HVREsOXs 2HVR» 5X_ SHOEV ISX 
2 DHSUMs BXs Z2NS0s 9X2 SHR Eo a 5X2 ZHSUMZ OX, 2HSD ) 
ICT=0 
ICTE 2d 
Iotos0 
SUMLD#0.0 
SUMUSQ20,0 
SUMRE=0,0 
SUMRES#0,0 
ONE ALANK CARD USED T2 SEPARATE DATA FUR DIFFERENT TARGETS 
4 CARD WITH ENC PUNCHED IN THE FIRST THREE COLUMNS 
WILL TERMINATE THE PROGRAM 
150 READ(5s10) IGNs TRHOs TSHNs THICK» ANGs PMASS » POT As VSPs VRESRMASSS HOTA 
TFC TRHO.LE.0.QIGSTO 200 
COMPUTE PROJECTILE CRISS=SECTIONAL AREA 
ARE.stP1¥(POTA/2,0)*4#2 
THE UNIT OF VELOCyYTY TN BE USED IS CM/SEC 
VS"VSP*100.0 
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amon 


aan 


160 


wae 


VR30.0 
DELV#10000.0 
THICKPoTHICK+DEL TAX 
THICKM2THICK=DEL TAX 
COMPUTE COEFFICIENTS AND OTHER QUANTITIES 
(NOTE + 9,867 CUNVERTS THE BRINELL HARDNESS NUMBER 
FROM KG/MM¥*#2 TQ DYNE/CMS#2 ) 


‘CA#9, BE7T#TBHNOC] 


CBsSQRT(9¢ BE 7THTBHN*TRHO) 4C2 
CC2TRHO¥C3 
QX84 OC 1L9C3-CQe42 
Q2QX*9, BE7T#TBHNZTRHO 
QOsSQRT(Q) 
QLECA+CAMVS eC CHV Sae2 
Q332,0*CC#VS§ 
Q52(COS(ANG*OTOR) )*#1.05 
Q9"1),5/CC#( PMASS/AREA) 
XTO=Q9%(ALOG(Q1L/C4)42-.64#CB/Q0*(ATAN(CB/Q0) eATAN( (Q34CB)/Q0))) 
XTO=XTO#QS 
WHEN XTO < THICK » THEN PENETRATION IS \NCOMPLETE «= VRa0.0 
LFAXTO LT eTHICKIGOTO 165 
“OMPUTE A FIRST ESTIMATE FOR THE RESIDUAL VELOCITY 
RATEN=ABS ((XTO@THICK)/THICK) 
LECRATIO GT LO) RATIO#ABS ( (XTO@THICK)/XTO) 
VRERATIO*VS 
TFLGOM=O 
TFLGP SO 
TREP 30 
TREPSTREP+) 
IFC IREP.GT.1L00)GUTU 910 
IF CIFLGP.EQ,1L.AND.TFLGM,2Q.1)G0TO 180 
COMPUTE THE PREDICTE) TARGET PLATE THICKNESS 
Q220 AF+CB¥VROCCHYVROR2 
Q4Gma7 .UeCCRVR 
XT2W9%(ALOG(Q1/Q2)42,9%CB/Q0%( ATAN( (044CB)/Q0)=ATAN( (Q346B8)/20))) 
ACCOUNT FOR OBLIQUE ANGLE IMPACT 
XT2XT¥*Q5 
([FCKXT.GT«THICKPIGOTO 170 
TF(XT.LTSTHICKM)GOTS 175 
FALLING THROUGH THE ABOVE IF STATEMENTS MEANS THAT TOLERANCE 
JN THE PLATE THICKNESS HAS BEEN MET WITH THE CURRENT VALUE 
FOR THE RESTOUAL VELOCITY 
GOTO 185 
VREVR#DELV 
IFLGP=#] MEANS YHAT THE VALUE FOR VR IS TOO LOW 
IFLGP #1 
GTO 160 
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175 
¢ 


185 


186 
187 


189 


195 


VReVR-CELV 
IFLGM=1 MEANS THAT THE VALUE FOR VR JS TOO HIGH 
1FLGM21 
GOTO 160 
DELV®DELV/2,0 
THE VALUE FOR VR HAS OVERSHOT THE TOLERANCE LeveL FOR 
THICKNESS MEANING THAT THE VALUE FOR DELV IS TOO LARGE 


“TFLGP #0 


IFLGM=0 
IF(DELV.LT.O-1)9G0TO 910 
GOTO 160 
VRAVR/10000 
COMPUTE DEVIANT AND RELATIVE FAROR AND CORRESPONDING SUMMATION 
DaVReVRE 
ICTORtCTOF) 
SUMD=SUMD+D 
SUMDSQsSUMNS§Q+Dea2 
RELFERR#1000,0 
IF(VRE.LE+6,0)GOTO 1a6 
RELERR#D/VRE 
GOTU 187 
IF CVReLE2O,0)RELERR#0.0 
IF (RELERReGE.500.0)G0TO 189 
ICTERICTE?)] 
CONTINUE © 
SUMRE =SUMRE*RELERR 
SUMRES@SUMRES+RELERRS*2 
WRITE (69195) IDNs TRHO, TRHN» THICK» ANG» PMASS es PDITAs VSP» RMASSSHN IAs 
VRE» VRafs SUMO, SUMDSQ, RELERR» SUMREs SUMRES 
FORMAT( LA sAGeFO,29 FoelsF8e3eF5,. LrFPO ede be3dsFB,leF 7,22 
FV .Se4FBe}eFlOCL, SFA. ) 
ICTHICT+1 
GAOTN 150 
CONTINUE 
, FIND THE MEAN, VARIANCE ANI) STANDARD DEVIATION OF 
“DEVIANTS ANO RELATIVE ERROR 
CTsiCTO 
CT i's IcThe1 
NBAR=SUMD/CT 
NVAR s(SUMDSQ=NDBAR#*F2eCT)/CTI 
DVSD=SQRT(OVAR) 
CTsIcTe 
CTLFICTE=1 
EBAR#SUNMRE/CT 
EVAR=(SUMRES~EBARS*2aCT)/CTL 
ESD-SQRTCEVAR) 
WRITE (62220) 1CTD 
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WRITE (6s220) 0BARSCVARSDSD 
WRITE (6s 220) ICTE 
WRITE (69215) EBAR,EVARZESD 


210 FORMAT CLHOs10xXs37HMEANs VARIANCE AND STANDARD DEVIATION 2X, 


1 


14HOF DEVIANTS = » 3F10.1 ) 


215 FORMAT(1HO,10X237HMEAN, VARIANCE AND STANDARD DEVIATION s2X> 


1 


* 

TUPAST 
1.0 
209 


17HOF RELATIVE ERROR » 3F10,5) 
PRINT NUMBER GF POINTS 
WRITE (69220) ICT 
FURMAT(1HOs10Xs ). BHNUMBER OF POINTS = » 15} 
IFC TONeEQe 6HEND 1GQTO 900 
GoTo 100 
WRITE (62905) 
FORMAT(1HOs30Xs LOHEND OF RUN ) 
sTop 
WRITE(6s915) ION, TREPsDELY 
FURMAT(1HOs L0Xs Abs 2X, 19HNID NOT CONVERGE IN » 13592X> 
LOHITERATIONS25X26HOELV = » F10.5 } 


GNTH 150 

END 

DATA 

UM ALLOY 
4.48 190.0 2127 -9 1.95 4759 567.84 821.21 
4.48 190.0 127 00 1.695) «6 759 1461.82 1294.18 
4.4¢° 190.0 318 -O 1.95 4759 880,26 390.40 
4.48 190.0 2318 20 1.95 .759 1355.45 1083.56 
404A 190.0 2635 -O 1695 6759 1491.08 683.36 
4.48 190.0 2635 .0 1.95 4759 1986.38 1127.15 
4.48 190.0 1.270 -2 1695 6759 2371.65 361.6) 
4.44 190.0 2127 ~)? 3.89 16013 798.85 672.69 


4.48 190.0 2127 19 3.89 16013 1932.66 957.58 
4.48 190,90 2318 -9 3.89 16013 620,27 500.79 
4.4 190.0 .318 .9 3689 16013 773.28 582,76 
6.4e 190.0 w319 -O 3,89 1.013 1499,01 1251.81 
4.4¢ 190.0 635 .) 3.89 16013 1505.71 774.50 
4.46 190.0 2635 .O 3.89 16013 1526.13 979.32 


4.44 190.0 2935 ~9 3,89 16013 2455.16 1367.33 
4.48 190.0 1.279 29 3.89 1.013 2551.79 1165.86 
4.4% 190.0 127 29 7.78 12267 640,99 563.44 
4.44 190.0 127 -9 7.78 16267 959,21 8746.78 
4.48 190.0 2318 .9 7.78 16267 Y7D,97 785.47 
4.46 190.0 2635 29) = 7.78 12267 149864,38 994.87 


7.78 135.9 24H 22 1.95 759 888,49 848.87 
7.76 135.0 2192 0 1.95 6759 1211.58 1914.98 
7.78 300.0 2318 9 1.695 759 1521.26 1196.34 
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7.7? 
10299) 


-10.9° 
-10.9" 
“10,90 


1869 


197 
198 


209 
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eoovo0o000nro 9020 200 


be de coe cel oe ell eel eel oll 


7.78 
7.78 
7.78 
7.78 
7.76 
7.78 
7.78 
7.76 
7.78 
7.7€ 
7.78 
7.76 
7.78 
7.78 


305.0 
393.0 
135.0 


135.0 © 


135.0 
300.0 
300.0 
300.0 
300.0 
300.0 
300.0 
300.0 
305.0 
332.0 


1.270 


1.93 
1.95 
3.89 
3.89 
3.89 
3.89 
3.09 
3,89 
7,78 
15.36 
15.56 
15.56 
15.56 
15.56 
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e759 

0789 
1.013 
3.013 
1.013 
12013 
12013 
1-013 
1026? 
1249) 
1249) 
10491 
1649) 
1049) 


1394.46 
609,90 
302,06 
393,50 
883.62 
879.65 

1466,09 

1466,09 
909,52 
916,53 

1425.53 

1432,56 

1556.00 

1660,55 


460.25 
967.89 
277.37 
236.03 
342.54 
583.69 
1164.36 
687.32 
690.37 
734.38 
1179.58 
1037.84 
1109.47 
759.56 


~10.00 
10,06 
10,0) 
"16-09 
“10.90 


“710.00 


<10.69 
#10,00 
-10,00 
#16,00 
710.09 
-10,60 
-10.00 
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PROGRAM USING THE THOR EQUATION TO PREVICT RESTOUAL VELOCITY 
GLUSSARY OF VARIABLES 


IDENTIFIES REQUIRED INPUT DATA 
IDENTIFIES INPUT DATA WHICH IS NOT REQUIRED 


ANG ~- THE ANGLE OF THE TARGET PLATE WITH RESPECT TG 
LINE OF FLIGHT = GBLIQUITY (DEGREES) 
eRerA = THE PROJECTED CROSS-SECTIONAL AREA OF PROVECTILE 
ON  MPaAcT 

n - THE DEVIANT 7% COMPUTED VALUE MINUS EXPERIMENTAL VALUE 
VBFR ~- THE AVERAGE (MEAN) VALUE UF THE DEVIANTS 
nS = THE STANDARD DEVIATION OF THE DEVIANTS 

QTOR ~ CONVERSION FACTOR ~ DEGREES TO RADIANS 
NVAR © THE VARIANCE OF THE DEVIANTS 

€1 = CUNSTANT BASED ON LEAST SQUARE FIT TO DATA 

12 ~ CONSTANT BASED ON LEAST SQUARE FIT TO DATA 

73 - COMSTANT BASED ON LEAST SQUARE FIT FC DATA 

4 ‘= CONSTANT BASED ON LEAST SQUARE FIT TO DATA 

5 - CONSTANT BASED ON LEAST SQUARE FIT TO DATA 

cBaAR = "HE AVERAGE (MEAN) VALUE OF THE RELATIVE ERROR 

ESt ~ THE STANDARD vEVIATION OF THE RELATIVE ERROR 

2 VAR ~ THE VARIANCE OF THE RELATIVE ERROR 
“DIA «© THE DIAMETER SF THE HOLE MADE IN THE TARGET (CM) 
en - INDEX TO COUNT NUMBER OF DATA CAROS FOR ONE TARGET 
ct) - INDEX COUNTER ON NUMBER OF POINTS FOR DEVIANTS 
"CTE ~ INDEX COUNTER ON NUMBER OF POINTS FOR RELATIVE FRROR 
DN - AN TOENTTFICATION NUMBER JR SYKBOL=*DESIGNATES SHIT NR, 
PDIA ~ DTAMETER OF THE PROJECTILE (CM) 
PMASS ~- MASS OF PROJECTILE (SRAMS) 
RELERR = THE RELATIVE ERROR DF COMPUTEO VS, EXPERIMENTAL 
RMOSS = THE RECOVERED PROJECTILE MASS (GRAMS) 

SUMN - THE SUM OF THE DEVIANTS 

SUMNSQ = THE SUM OF THE DEVIANT SQUARED 

SUMRE — THE SUM OF THE RELATIVE ERRUR 

SUMRES «- THE SUM OF THE RELATIVE ERROR SQUARED 

Y Bhi ~ THE TARGET BRINELL HARONESS NUMBER (KG/SQ MM) 
THICK -— THE THICKNESS QF THE TARGET PLATE (CM) 

TRHI) = THE nENSYTY OF THE TARGET PLATE (G/CC) 

VR = PREDICTED RESIDUAL VELOCITY (M/S) 

VRE = THE EXPERIMENTAL RESIDUAL VELOCITY (M/S) 
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ae VSP - THE EXPERIMENTAL STRIKING VELOCITY (M/S) 


VATA PISDTOR / 34141892634» 060174532925 / 
1G FORMATC Abs 2Fbels FS .3eF5.le2FO, 32 2FOele2FT,2) 
11 FORMAT(5F10,3) 
12 FORMAI(10A6) 
20 FORMATCLHOsLOXs 1 3HEXPONENTS ARE » 5F1O.3 / ) 
100 ‘CONTINUE 
WRITE( 63103) 
194 FOR AAT (IHL) 
WRITE (69106) 
196 FURAATCLHOs20XeG2HSAMPLE SJUTPUT FOR THOR EQUATION PREDICTING » 
L 18H RESIDUAL vELOcITy / ) 
READ CARD WHICH IDENTIFIES TARGET MATERIAL 
REAU(CSsL2Z9TGT Ls TGT2sTGT3 
WRITE(6912)°GT1L,TGT2,TGT3 
REAI(CSeLLIELZ EZ, E32 EGSVES 
WRITE(6220)E1,E2, 63, E4, ES 
WRITE( 69110) 
119 FOR 1AT(1HO¥4Xy39HNR TRHO TBHN THICK ANG MASS PDIA 246X> 
1 ?HVS92%s L2HRMASS HDIA »SXe3HVRE2OX>2NVR»e 5SXs SHOEV IS Xe 
? ZHSUMs BXs 2HS Qs 4X9 4HR eG 0d 5X9 SHSUM9OXs ZHSQ ) 
rT. 20 
oTisG 
SUMVEO LO 
SUM2§U20.0 
SUMRE=0,0 
SUMRES=0,0 
1GT26 
INE BLANK CARD VSEN TO SEPARATE DAIA FUR DIFFERENT TARGETS 
» GARG WITH &ND PUNCHED IN THE FIRST THREE COLUMNS 
WILL TERMINATE THE PROGRAM 
15Q REA 5s LOY TON, TRHO2 TAHNs THICK» ANGs PMASS » PUTA, VSP, VREsRMASSSHIIA 
fF¢TRHO LE,G.0)GOTU 200 
VOMPUTE PROJECTILE CROSS-SECTIONAL AREA 
AREAsPT¥(POITA/2.0)%%2 
THF UNIT OF VELOCITY TO BE USED IS CM/3EC 
VS 3VSP*100,0 
91210,0%8E) 
Q2=( THICK*AREA) ee b2 
Q3J="MASSHHEY 
Q42(1.,0/C0S8 (ANGE TOR) ) HES 
Qoev§xees 
-GMPUTE THE PREDICTE! RESIDUAL VELUCITY 
VRStVS=Q1402%Q3404%05)/190,0 
-GMPUTE OEVIANT aNd 2ELATIVE ERROR AND CORRESPONDING SUMMATION 
WeVReVRE 
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SUMO SSUMO4D 
SUMUSQ=SUMNSQ4D%82 
1CTHelCTO+) 
RELERR#1000,0° 
TF(VRE.LE*O-0)G0TD 1846 
RELERR=D/VRE 
GOTO 187 

Lib TECVR LEO. Q)RELERR «0.0 

187 LFCRELERR.GE.800.0}G0TO 189 
ICTHSICTE+1 

189 CONTINUE 
SUMRE* SUMRE*RELERR 
SUMRES*SUMRES+RELERRW#2 


WRITE(69 195) LDN» TRHG, TBHN» THICK, ANGs PMASS + POL As VSP>RHASS sHD1 As 


VRE» VRe Ds SUMDs SUMDSQ:; RELERRs SUMREs SUMRES 
LUD FORMATCLH sAOsFO.2eFbeLr FO SoFS,LaFO ater Oude FO. lor lols 
1 FT. Ba4FOolsFlOels3FA.3 ) 
ICT ICT+} 
GOTO 150 


200 CONTINUE 


FIND THE MEAN, VARYANCE AND STANDARD DEVIATION OF 
DEVIANTS AND RELATIVE ERROR 
cTsicto 
CTLsICTN|) 
OXARSSUMD/CT 
JVAK®(SUMDSQ=DBAR*2eCT)/TT) 
OSN=SGRTEDVAR) 
cYsICTE 
CTLHIcTeE=) 
FNXAR=ZSUMRE/CT 
EVAR = (SUMRES<EBARS*2aCT)/CTL 
ESD«SQRTCEVAR) 
WRITE (G,22G)ICTO 
WRITE CGs210) DBARSNVARZDSD 
WKITE CO, 220) CTE 
WRITECOs 215) EBAR,EVARSESO 
210 FURLALCTHOs 10X23 7HMEAN» VARIANCE AND STANDARD DEVIATION »2X> 
‘ L4HUF DEVIANTS © » AF10.1 ) 
215 FNRMATCLHOs LOXe37HMEAN, VARIANCE AND STANDARO DEVIATION 52k, 
1. U7HQE RELATIVE ERROR » 3F 10.5) 
PRINT NUMBER OF POINTS 
WRITE 69220) 1CT 
220 FURMAT(CLHOs,1LOXs 1 BHNUMBER OF POINTS = » 15) 
IF (SDt.FQ. #HEND 1GOTN 900 
GUTE L0G 
9OO WRITE (6,905) 
9G5 FOR AT(1HOs30X2 1. OHEND OF RUN ) 
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END 243 

* DATA 444 
TTTANTUM ALLOY ’ . 

4,885 1.103 -1,095 1+369 0.167 

1.0 4.48 190.0 127 0 1.99 73% $67.84 $21.21 1.88 

260 4,48 190.0 e127 -O 1.93 759 1461.82 1294.18 091 

36% 4648 190.0 2318) =,0) 1.98 «6759 «880.26 390,40 -10.00 

4.9 4.44 190.0 2316 0 1.95 2759 1355.45 1683.56 1.63 

be 4.468 190.0 635 0 3.95 6759 1491,08 683.36 1.73 

oO 4.48 390.0 eo (§ .0O 1.935 759 1966.38 1127.15 56 

Te 4648 190.0 1.270 -O 1.95 759 2371,65 381.61 ~10.00 

6,0 6.46 190.06 2147 .6 3.89 1.013 798,68 072.69 3.83 


Pe) 4.48 1390.0 2127 Pa?) 3.89 1.013 1032,66 957.68 1.72 
10.9 4.49 190,6 318 -O 3.89 16013 620,27 500.79 3.83. 
Yye0 4.45 190.0 318 .O 3.869 1.013 773.28 582.78 3.81 
12.0 4.48 196.9 38 -9 «3.89 1.013 1499.01 1251.61 2243 
13.0 4.468 190.9 638 20 3,89 1.013 1505.71 774.50 3624 
14.f 4.48 190,06 2635 2D 3,69 1.013 1526.13 979,32 ~10.00 
15.0 4.42 190.0 2635 oO 3,89 1.013 2455.16 1367.33 el2 
16.0 4.46 190.9 ~270 .O 3.89 1.013 2551.79 1168.86 057 
7. = 4 19020 2127 .O 7.78 16267 640,99 561,44 7.72 
1H.eO «64488 190.9 127 oO «7.78 16267 959,21 874.78 7.72 
19,0 4.48 190.0 2348 0 = 7.78 Le267 975.97 785.47 772 
20-9 4.48 190,06 7695 0 7.78 1.267 1484.36 994.87 =30,00) 


~~ 


§,699 0.649 @V 2945 0.989 0.619 

1.9 7.78 135.9 2046 .o 1.95 759 888,49 848.87 ~10,.00 
2% 7.78 135.0 2192 20 1,98 759 1211,58 1014.98 ~10.00 
369 97.7% 300.6 2318 20 «1698 6759 1521.26 1196634 ~10.00 
4% Fath 303.6 695 20 1495 6759 1394.46 460.28 ~10.90 
5.0 7,78 393,0 2152 29 1.495 6759 609,90 367.89 ~10.90 
+o) 7.7b 335.0 2946 ~O 3.89 16013 302,06 277.37 -10.9n 
7.0 7.78 £35.0 .1°2 oO 3,89 34013 393.50 256.03 ~10.00 
Hel eT» 135.9 2316 .O 3.49 146013 883.62 4342.54 ~10.00 
Oo F,78 306.0 o3h8 ~? 3.89 1.6013 879,65 $83.69 ~10.00 
Oo t.76 300,90 23h8 .O «3.89 1.6013 1466.09 1164.34 -:0.00 
Oo 4.7» 300.0 »638 oO «3,89 1.013 1464.09 687.32 ~10.00 
Oo 1.72 300,0 wala 7 «#87278 16267 909,52 690.37 ~10.00 
f 7.7 2080.9 2318 -O 13.56 16491 916,53 754.38 ~10.90 
144.0 7.7 300.0 og18 : 
tS. P.78 360.0 2635 : 
Jb 7.78 3205.0 2635 ® 
17,0 7,78 332,39 1.278 : 


i) 

O 45.56 16491 1425.55 1179.58 -10.00 
@ 15.56 16491 1432.56 1037.86 -10.00 
0 15,56 16491 1556.00 1109.47 ~10.00 
0 16.56 16491 1660.55 789.56 -10.0¢ 
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PROGRAM USING THE 2/F EQUATION TO PREDICT PLATE THICKNESS 
SLOSSARY OF VARIABLES 


TOENTIFIES REQUIRED INPUT DATA 
IDENTIFIES INPUT DATA WHICH JS NOP REQUIRED 


NG ~ THE ANGLE GF THE TARGET PLATE WITH RESPECT TO 
LINE OF FLIGHT = QOBLIQUITY (DEGREES) 

AREA = THE PROJECTED CRUSSSECTIUNAL AREA OF PROVECTILE 
UN IMPACT 


oy ~- CUNSTANT BASED ON LEAST SQUARE FIT TO THOR DATA 

Le - CONSTANT BASED ON LEAST SQUARE FIT TO THOR DATA 

“3 ~ CONSTANT BASED ON LEAST SQUARE FIT TO THOR DATA 

iy) = THE DEVIANT «= COMPUTED VALUE MINUS EXPERIMENTAL VALUE 
RAR ~ THE sVERAGE (MEAN) VALUE OF THE DEVIANTS 

usp =~ THE STANDARD DEVIATION OF THE DEVIANTS 

OTUR = CONVERSION FACTOR = OEGREES TO RADIANS 

VaR “= THE VARTANCE 2F THE DEVIANTS 

BAR - THE AVERAGE (MEAN) VALUE UF THE RELATIVE ERROR 

“$§L ~ THE STANDARD DEVIATION OF THE RELATIVE ERROR 

eVrkK = THE VARTANCE UF THE RELATIVE ERROR 

et1A » THE DIAMETER OF THE HOLE MADE IN THE TARGET (CM) 
"eT ~ INDEX TO COUNT NUMBER OF UATA CARDS FOR ANE TARGET 
ope ~ AN TDENTIFICATION NUMBER UR SYMBOL@DESIGNATES SHUT NR, 
"Tr 1a ~ DIAMETER QF THE PROJECTILE (&M) 

ruLSS «— MASS OF PROJECTILE (GRAMS) 

oEL ERR = THE RELATIVE ERROR OF COMPUTED VS, EXPERIMENTAL 
“MiS§ = THE RECOVERED PROVECTILE MASS (GRAMS ) 

el) ~ THE Sub CF THE DEVIANTS 

SUrDSO = THE SUM OF THE DEVIANT SQUARED 

SUMRE = THE SUM MF THE RELATIVE ERROR 

SUMRES ~ THE SUM MF THE RELATIVE ERROR SQUARED 

TREN ~ THE TARGET ARINELL HARDNESS NUMBER (KG/SQ MM) 
THICK «© THE THICKNESS OF THE TARGET PLATE (CM) 

TRHO ~ THE DEASITY OF THE TARGET PLATE (G/CC) 

vet - THE EXPERIMFNIAL RESIQUAL VELUCITY (M/S) 

VSF ~ THE FKPERIMENTAL STRIKING VELOCITY (M/S) 

xT ~ THE PREDYCTED TARGET PLATE THICKNESS (CM) 


WAT’ Pla DTOR / 3,1415892¢54s 960174532925 / 
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Seu utwae 


Ana 


10 FORMAT( AOs2FOelsFO.32F3 eis 2FOe3s2F8.i92F 7.2) 
12 FORMAT(10A6) 
C1*«0,70 
€220,.23 
€38u,50 
100 CONTINUE 
WRITE (69103) 
103 FORMAT( 1H) 
WRITE (63106) 
106 FORMAT(LHOs20XeS2HSAMPLE OUTPUT FOR Z2/F EQUATION PREQICTING » 
1 16h PLATE THICKNESS / ) ; 
READ CARD WHICH YOENTIFIES TARGET MATERIAL 
READ (5912) 76T1LsIGT2s 7G73 
WRITE(6212)7GT1L,TGT2,TGT3 
WRITE(6,110) 
110 FORMAT(}HOs4Xs31HNR TRHO TBHN ANG MASS POIA » 6Xe2HYSs6Xe 
) 24HVR RMASS HOTA THICK » 6Xe2KXT95Xs SHDEVs SX0 
2 BHSUMs Xs 2H8Qs 4X5 SHR OE 2 5X2 9HSUMZOXs 2HSQ ) 
SUM 20.0 
SUMDS$Q80,.0 
SUMRE30,0 
SUMRES=0.0 
ICT#=0 
ONE BLANK CARD USBO TO SEPARATE DATA FOR DIFRGRENT TARGETS 
A CARD WITH END PUNCHED IN THE FIRST THREE COLUMNS 
WILL TERMINATE THE PROGRAM 
150 READ (5510) LON» TRHQs TOHNs THICKs ANGs PMASSs POL As VSPaVRESRMASSSHOTA 
IF (TRHOLEC.0)GOTO 200 
COMPUTE PROJECTILE CROSS@SECTIONAL AREA 
AREA=P16(PDJA/2.0) 482 | 
THE UNIT OF VELOCITY TO BE USED I§ CM/SEC 
VSsVSP*100.6 
VReVRE*100.0 
COMPUTE COEFFICIENTS ANDO OTHER QUANTITIES 
(NOYE + ©@,8E7 CONVERTS THE BRINELL HARONESS NYMeaR 
PRON KG/MM#82 TQ DYNE/CHee2 ) 
CA=9. BE 7eTBHNGC] 
CBuSORT (Oe OS 7STBHN®TRHO) C2 
CCeTRHOSCS 
QX24.08C 166 3=C2ae2 
Q2QX69, SETOTBHNOETRHG 
Q0sSQRT(Q) 
QLeCAsCBeVvEeCceyses2 
Q2sCA+CBOVASCCeVReS2 
Q3=2.,07CCOV§ 
Q452,08CCOVR 
Q52(COS( ANGSOTOR) ) #83 .05 
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Q980,5/CC®(DMASS/ AREA) 


Cc COMPUTE THE PREUICTEN TARGET PLATE THICKNESS 
XT2Q9%(AL0G(Q1/G2)4*2,96CB8/Q0«( ATAN( (Q49CB) /Q0) -ATAM( (034068) /00))) 

C “CCOUNT FUR OQaLiQue ANGLE IMPACT 
XTsxT#Q5 

C COMPUTE OEViANT AND RELATIVE ERROR AND CORRESPONDING SUMMATION 
NeXT=THICK 


“SUNG #SUMO4D 
SUM )SQ*SUMDSQ+D422 
RFLERR=D/THICK 
SUMRE=SUMRE@RELERR 
SUMRES SSUMRES+REL ERR&*2 
WRITE (02155) 10Ny TRAD, TBHN2 ANGs PMASS» POT AS VSP» VRESRMASS* HDI A> 
1 THICKsXTeD 2 SUMD» SUMNSQRELERRs SUMRE »SUMRES 
LGTalc Tel 
GUT 150 
LoS BRIRMATOCLH » Aba Fo. 22F6.1525.19F6.29F6,392F9,152F7,2,8F8, 3) 
2oo cYsict 
oTpe ie Te) 
+ IND THE MEAN, VARSANCE AND STANDARD DEVIATION OF 
VEVIANTS AND RELATIVE ERRUR 
JV ARRSUMD/CT 
VAR ® (SUMOSQ=<NBAR®*28CT) /CT) 
DSO-SARTCOVAR) 
- JAR BSUMRE/CT 
KVAK&(SUMRES“EBAR**2maCT)/CT1 
“50>. SQRTCEVAR) 
ARLE (62210) 0BARSDVARINSD 
WRITE COs 215) EBARZEVARZESD 
217 FE IRAATCLHOs, LORS 3THMEAN, VARIANCE AND STANOARD DEVIATIUN 22k, 
1 14H0F OEVIANTS ® » 3F10.5 } 
21% FYURUAT(IHOs LOXs37HMEAN, VARIANCE AND STANUVARO DEVIATION »2X> 
\ ‘7RUF RELATIVE ERROR » 3F10.5) 
4 SRINT NUMBER OF POINTS 
WRITE(Os220)1CT 
220 FOR ATCIHOs1L0Xs LAHNUMBER JF POINTS = » 15) 
TFC IDN.LEQs OHEND 160TH 900 
GUT LOO 
Ine WKTTE (00905) 
YNS FOR. ATCLHOs30X, 1OHEND OF RUN ) 
ST? 
FO 
* waTlA 
TIT Gtg4* ALL AY 
1e% 44H 190.0 2)27 -O 1.95 4.759 567.84 521.21 1.88 
20% 4.4" 190920 2127 = ,9 1,95 0759 1464-82 1296.18 09) 
1. 4.408 190.0 2318 -2 1.95 759 580,26 3890640 -19.9°% 
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oo (i Sees 
a ta ie 
ug SMU Bas: 


ed 
Nw CODA RMU-AS 
e e e oe . e e ee @ 
9oq 0090 9O0090090 


RHA 


eecees8 ee ee seee 
ome I> Ihe Ibe Ie Ee | 


oODUDII SD 


14. 


4.48 
4.48 
4.48 
4.48 
4.48 
4.48 
4.48 
4.48 
4.48 
4.48 
4e48 
4.48 
4.48 
4.48 
4.48 
4.48 
4.48 


7.78 
7.78 
7678 
7.78 
7.78 
7.78 
7.78 
7.78 
7.78 
7.7% 
7.78 
7.78 
7.78 
7.78 
7.78 
7.78 
7.78 


318 .0-- 1695 
435 ~? 1.95 
2695 2 1.95 
1.270 Oo «1.95 
.L27 .o 3,89 
127 .0 3,89 
318 -O 3,89 
2314 9 #%3,89 
.2l8 .o %3.,89 
635 .o 3,89 
2035 -o 3,89 
2635 Oo 3,89 
1.279 Oo 3.89 
~l27 Oo 7,76 
127 Oo 7.78 
318 sO R78 
2635 .0 7,78 
246 .o «1,95 
2152 Oo 1,95 
2318 Oo 1,95 
2935 oO 1,95 
2152 Oo 1,95 
2046 .2 3,89 
2152 0 3,89 
2318 .o 3,89 
2318 .o 3,89 
31a .O 3,89 
2635 Oo 3,89 
2318 .O 7,78 
318 -0 15,56 
318 .o 15,56 
~635 .o 15,56 
.635 20 15,56 
1.270 .2 15,56 
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2759 

0759 

0759 

0759 
1.013 
1.013 
1.013 
1.013 
1.013 
12013 
1.013 
1.013 
12013 
1.6267 
1¢267 
16267 
1.267 


0759 

2759 

0759 

0759 

0759 
16013 
12013 
12018 
16013 
1.013 
1-013 
1.267 
1249) 
10491 
1.49) 
1449} 
1649} 


1355.45 
1491.08 
1986,38 
2371.65 
798,88 
1032.66 
620.27 
773.28 
1499,01 
1505,71 
1526.13 
2455.16 
2551,79 
640,99 
959.21 
975,97 
1484,38 


888.49 
1211,58 
1521,26 
1394.46 

609.90 

302,06 

394,50 

883,62 

879,65 
1466.09 
1466,09 

909,52 

916,53 
1425,55 
1432,56 
1556,006 
1660,55 


1083.56 
683.36 
1127.15 
381.61 
$72.69 
957.68 
500.79 
382.78 
1251.81 
774.50 


979.32 : 


1367.33 
1165.86 
5361.44 
674.78 
785.47 
994.87 


848.87 
1014-98 
1196.34 

460.25 

367.89 

277.37 

256.03 

542.54 

583.69 
1164.34 

667.32 

690.37 

754,38 
1179.58 
1037.84 
1109.47 

759.56 
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-J0.90 
-10299 
“10.99 
710690 
~1009 
#10.9> 
-10e9) 
-10.9) 
710.99 
710-9 
“10.90 
-10.99 
#10299 
~10.0% 
-10.95 
“10.99 
“10.99 
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rROGRAM USING THE THOR EQUATION TQ PREDICT PLATE THICKNESS 


SLUSSARY OF VARIABLES 


IDENTIFIES REQUIRED INPUT DATA 
IOBNTIFIES INPUT DATA WHICH IS NOT REQUIRED 


\NG - THE aNGLE QF THE TARGET PLATE WITH RESPECT TO 
LING OF FLIGHT = OBLIQUITY (OEGREES) 
ReA = THE PROJECTEO CROSS—SECTIONAL AREA DF PROJECTILE | 
ON yMPact 
THE DEVIANT = COMPUTED VALUE MINUS EXPERIMENTAL VALVE 
THE av’ tRAGE (MEAN) VALUE UF THE DEVIANTS 


$2 ~ THE © ANODARN DEVIATION OF THE DEVIANTS 

JTAR = CON STON FACTOR = OEGREES TO RADIANS 

YVAR © THE TANCE OF THE DEVIANTS 

=y - CONS, .NT BASED CN LEAST SQUARE FIT TO DATA 

£2 - CONSTANT BASED OV LEAST SQUARE FIT TO DATA 

a) - CONSTANT BASED On LEAST SQUARE FIT TO DATA 

m4 ~ CONSTANT BASED ON LEAST SQUARE FIT TO DATA 

bs - CONSTANT BASED ON LEAST SQUARE FIT TO DATA 

rBAR = THE AVERAGE (MEAN) VALUE OF THE RELATIVE ERROR 

1 $0 - THE STANDARD DEVIATION OF THE RELATIVE ERROR 

rVAR ~ THE VARIANCE OF THE RELATIVE ERROR 

4OrA ~ THE NIAMETER GF THE HOLE MADE IN THE TARGET (CMM) 
1cT - INDEX TO COUN! NUMBER OF DATA CAROS FOR ONE TARGET 
iON - AN INENTYFICATION NUMBER UR SYMBOL-DESIGNATES SHIJT NR, 
POIA = DJAMETER JF THE PROJECTILE (CM) 

VMASS = MASS MF PROYECTILE (GRAMS) 

SEL ERR = THE RELATIVE ERROR OF COMPUTED VS, EXPERIMENTAL 
KMASS © THE RECOVERED PROJECTILE MASS (GRAMS) 

SUMU ~ THE SUM QF THE DEVIANTS 

SUMDSQ = THE SUM OF THE DEVIANT SQUARED 

SUIRE = THE SUM QF THE RELATIVE ERROR 

SUMRES = THE SUM OF THE RELATIVE ERROR SQUARED 

TAHN = THE TARGET BRINELL HARDNESS NUMBER (KG/SQ MM) 
THICK = THE THICKNESS OF THE TARGET PLATE (CM) 

TRHO = = “THE pEMSyTy OF THE TARGET PLATE (6/C€¢) 

VRE ~ THE EXPERIMENTAL RESIDUAL VELOCITY (M/S) 

VSP - THE EXPERIMENTAL STRIKING VELOCITY (M/S) 

xT «= THE pRENYCTED TARGET PLATE THICKNESS (CM) 
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OATA PIsOTOR / 3413415992654 060174532925 / 
FORMAT ( AOs2FGole FO .32F 5. lv 2ZFOe de 2F Bele 2F7.2) 
FORMAT(S5F10,3) 
FORMAT(10A6) 
FORMATCIHOs1LOXe13HEXPONENTS ARE » 5F10.3 / ) 
CONTINUE 
WRITE( 60103) 
FORMAT(IHL) 
WRITE( 69106) 
FURMAT(1LHOs20X942HSAMPLE OUTPUT FOR THOR EQUATION PREDICTING » 
LOH PLATE THICKNESS / ) 
READ CARO WHICH JOENTIFIES TARGET MATERIAL 
REA Se L29TGTLSTGT2s,t6T3 
WRITE (6212) TGT1s,T7GT2,TGT3 
READ( Ds LLIELZE2, E32 EG2ES 
WRITEC Ss 20)E1,E2,E3,E4,ES 
WRITE (62110) 
FORMAT(1LHOs4X%s3)HNR TREO TRHN ANG MASS PDIA » 6Xa2HVS96Xs 
24HVR RMASS HDTA THICK » 6Xe2HXT25Xs3HDEV25Xe 
BHSUMS OXs 2HS Qs 4X _ SHR eE oa SX SHSUMZ OX 2HSQ ) 
SuMD320,0 
SUMS Q20,0 
SIMREz9.0 
SUARES20,0 
1eT=0 
"Nr BLANK CARD USED TO SEPARATE DATA FOR DIFFERENT TARGETS 
A CARD WITH END PUNCHED IN THE FIRST THREE COLUMNS 
WILL TERMINATE THE PRIGRAM 
REAM (5,10) LON» TRHUe TBHN»s THICK» ANGAPMASS PDT As VSP» VRE RMASS ®HOTA 
LF CTRHI.LE.O.QIGNTY 200 
/TIMPUTE PROJECTILE CRISS-SECTIONAL AREA 
AREAEP IRC PDIA/2,03 *¥2 
THE UNIT OF VELUCLTY TO BE USED IS CM/SEC 
v52VSP4100.0 
VR SVRE*1O0.0 
QL=lLO.0%*E1 
Qs=EPMASS#¥ED 
Q42(1.0/C05 (ANGHDTOR) ) ##E4 
Qo aVS#eES 
V.UMPUTE THE PREDICTED TARGET PLATE THICKNESS 
XT=((VS=VR)/(O18G3¥*Q48Q5) ) ¥*(1.0/E2)/AREA 
COMPUTE DEVIANT AND RELATIVE ERROR AND CORRESPONDING SUMMATION 
He XT=THICK 
SUMOeSUMOSE 
SUMISG=SUMDSQ40e#2 
RELERK=N/THICK 
SUM2E= SUMRE+RELERR 


1 


158 
290 


SUMRES*SUMRES+RELERR&S2 

WRITE (65155) LDNa TRHO, TBH» ANGs PMASS sPOT As VSP, VREsRMASSsHDIAs 
THICKI XT» Ds SUMDs SUMDSQ/RELERR» SUMREsSUMRES 

UCTHICT41 : 

GUTO 150 

FORAATIIH sAOs FO, 29 FG ele F 5.19 FG, 20F6,322F8,152F7,220F8, 3) 

CrslcT 


CTUtICT=1 


ocr. 
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cIND THE MEAN, VARIANCE AND STANDARI DEVIATION OF 
VEVIANTS AND RELATIVE ERRNR 

DNAR=SUMD/CT 

NVAKR®(SUMDSQ=DBAR#*2eCT) /CT1 

NSDV25QRT(DVAR) 

UEBAR#SUMRE/CT 

EVAR= (SUMRES-EBARe*2uCT) /CT] 

ESD=SQRT(EVAR) 

WRITE (69210) 0BARSOVARSDSD 

WRITE (Oe 215) EBARZEVARSZESN 

FURMAT(CLHOs LOXe397HMEAN, VARIANCE AND STANDARD DEVIATION »2X, 
14HOF DEVIANTS = » 3FILC.5 ) 

FORMATCLHOsL0Xs37HMEAN, VARIANCE AND STANDARD DEVIATION »2X, 
}7HOF RELATIVE ERROR » 3F10.5) 
PRINT NUMBER nF POINTS 

WRITE (69220) ICT 


» FORSATCINOs 10% *¢18HNUMBER JF POINTS = » 15) 


US CID EQe OHEND YGNTN 900 
SUT Ay 
WRITE (69905) 


FOR (ATC LUHOs30Xs LOHEND OF RUN ) 

stn: 

END 

ATA 

Meh SLLOY 

848 16103 “1.095 16369 0.167 


44 190.0 le? 2 1.95 759 567,84 52).21 1.98 
4.44% 190.0 127 2? 2.95 6759 1461.82 1294.18 oF) 
4.40 170.0 319 oO 1.95 6759 880.26 390.40 -10.90 
4.4" 190.0 232i o) 1.55 0759 1353.45 1083.56 1263 
44> 190.0 2635 2 1.95 6759 1491.08 683.36 1.73 
4.445 190.0 2935 9 = 1,95 2759 1966,3G 1127.15 93% 
4.45 190.0 1.2/9 oO 1.95 759 2371.65 381.61 -10.9) 
4.4F 190.y 127 Pa?) 3,89 1-013 798,83 672.69 32593 
464+ 190.0 aver Oo 3,89 146013 1032,66 957,68 1.7? 
4.4" 190.0 .3te .O 3.89 16013 520,27 500.79 3.83 
4.48 190.0 2318 ~2 3.89 1.013 773,28 582.76 3.81 
4.4: 190.6 318 Oo 3.89 1-013 1499,0) 1251.81 2.43 
4.4% 190.0 2535 .O 3,89 1.013 1565,71 774.50 32% 
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>2%2S2D 


190.0 
190.0 
190.0 
190.9 
190.90 


4 190.0 


190.0 


0.889 
135.0 
135.0 


4 300.0 


305.0 
393.0 
135.0 
135.0 
135.0 
300.0 
300.9 
300.0 


* 300.0 


300.0 


& 300.0 


300.0 
305.0 


» 332.0 


3,89 
3.869 
3.89 
7.78 
7.78 
7.78 
7.78 


1.013 
1-013 
1.018 
16267 
10267 
10267 
1.267 


0.989 


1,95 
1,95 
1.95 
1.95 
1.95 
3,89 
3.69 
3.89 
3,89 
3.89 
3.89 
7,78 
15.56 
15.56 
15.56 
15,56 
15.56 
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0759 
0759 
«759 
«759 
0759 


1.013 


1.013 
12613 
1eO18 
1.013 
1.013 
1.267 
12491 
1249) 
1491 
1491 
1.491 


1526.13 
2453.16 
2551.79 
640.99 
959.21 
975.97 
14864,36 


0,019 
868,49 
1211.56 
1521,26 
1394,46 
609,96 
302,06 
393,50 
A83,62 
879,65 
1466,09 
1466.09 
969.52 
916,53 
1425,55 
1432,56 
1556.90 
1660,55 


979.32 
1367.33 
1165.86 

$61.44 

874.78 

785.47 

994.67 


648.87 
1014-98 
1196.34 

460.25 

367.89 

277.37 

256.03 

$42.54 

$83.49 
1164.34 

687.32 

690.37 

754.38 
1179.58 
1037.84 
1109.47 

759.56 
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LIST OF SYMBOLS 
A projectile cross-sectional area at impact tom”) 
a Acceleration 
C projectile shape factor 
Ci pmpirical constants gaezis 2 and 3) 
F Resistive force 


u prinell hardness of target plate (aynes/em” 


t 

Ky A coefficient z Uy H, 

Ky A coefficient 2 Cc, Yo, HY 
K, A coefficient = C, O, 

m Mass 


Bi) projectile mass (grams) 


q A giscriminant ~ 4K Ky - 4 
~ Time to penetrate to depth x (sec) 
+ Time 
Velocity 


N 
V pesidual yetocity (cm/sec) 


ca] 


Vv Striking yelocity (cm/sec) 


w 


Ne Velocity at depth x 
Distance into the target measured from impact surface 
x, Target plate thickness (em) 


P, pensity of target piate (g/ce) 
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